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This  report  summarizes  the  progress  of  four  technical  investigations  conducted  during 
CY  90.  Although  reasonable  efforts  have  been  made  to  ensure  the  reliability  of  the 
data  presented,  it  must  be  emphasized  that  this  is  an  interim  progress  report  and  that 
further  experimentation  and  analysis  may  be  performed  before  the  conclusions  from 
any  of  these  investigations  are  formally  published.  It  is  therefore  possible  that  some 
of  the  observations  presented  in  this  report  will  be  modified,  expanded,  or  clarified 
by  our  subsequent  research^. 
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Abstract 

This  report  documents  the  technical  progress  in  the  four  investigations  which 
maAe  up  the  project  “Support  of  Research  Projects  for  Electrical  Energy 
Systems”,  Department  of  Energy  Task  Order  Number  137,  funded  by  the  U.S. 
Department  of  Energy  and  performed  by  the  Electricity  Division  of  the  National 
Institute  of  Standards  and  Technology  (NIST).  The  Erst  investigation  is 
concerned  with  the  measurement  of  magnetic  fields  in  support  of  epidemiogical 
and  in  vitro  studies  of  biological  field  effects.  NIST  cohosted  a workshop 
on  exposure  and  biological  parameters  that  should  be  considered  during  in 
vitro  studies  with  extremely  low  frequency  (ELF)  magnetic  and  electric  fields. 
Additionally,  equations  were  developed  to  predict  the  magnetic  field  in  a 
pajaUel  plate  magnetic  field  exposure  system.  An  IEEE  standard  prepared  at 
NIST  on  measuring  dc  electric  fields  and  ion  related  parameters  was  approved 
and  published.  Various  site  visits  were  made  to  characterize  the  electric  and 
magnetic  fields  in  biological  exposure  systems.  The  second  investigation  is 
concerned  with  two  different  activities:  the  detection  of  trace  levels  of  S2F10 
in  SFe  and  the  development  of  an  improved  stochastic  analyzer  for  pulsating 
phenomena  (SAPP).  The  detection  of  S2F10  in  the  presence  of  SFe  using  mass- 
spectrometric  detection  coupled  to  a gas  chromatograph  is  difficult  because 
of  the  similar  mass  spectra.  Enrichment  techniques,  capable  of  sub  parts 
per  billion  (ppb)  detection  in  the  case  where  the  background  gas  is  not  SFe, 
are  unsuitable  for  this  application.  A technique  is  described  that  enables  the 
detection  of  S2F10  in  gaseous  SFe  down  to  the  ppb  level  using  a modified  gas 
chromatograph-mass  spectrometer.  The  SAPP  has  been  improved  to  allow 
direct  real-time  measurements  of  conditional  pulse- amphtude  distributions  of 
higher  order  than  was  previously  possible.  The  new  system  was  applied 
to  an  investigation  of  the  stochastic  behavior  of  negative  corona  (Trichel 
pulses)  and  the  effect  of  a dielectric  barrier  on  these  discharges.  The  third 
investigation  is  concerned  with  breakdown  and  prebreakdown  phenomena  in 
liquid  dielectrics.  The  activity  reported  here  was  a study  of  negative  streamers 
preceding  electric  breakdown  in  hexanes.  Using  the  image  preserving  optical 
delay,  the  growth  of  the  streamers  associated  with  partial  discharges  at  a 
point  cathode  are  photographed  at  high  magnification.  Simultaneous  discharge 
current  measurements  enables  a detailed  description  of  the  temporal  and  spatial 
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development  of  streamers  and  provide  a basis  for  the  evaluation  of  models  for 
the  initiation  of  negative  streamers.  The  last  investigation  is  concerned  with 
the  evaluation  and  improvement  of  methods  for  measuring  fast  transients  in 
electrical  power  systems  such  as  might  be  associated  with  an  electromagnetic 
impulse.  The  new  draft  of  IEC-60  recommends  the  use  of  an  independent 
reference  measurement  system  to  verify  the  performance  of  the  impulse  voltage 
measuring  system  under  test.  A compact  resistive  divider,  NIST4,  was  designed 
for  this  purpose.  It  is  anticipated  that  this  divider  together  with  some  Kerr 
electro-optical  devices  will  be  used  as  the  reference  system  at  NIST.  The  design 
details  of  NIST4  and  its  measurement  capabilities  are  presented. 
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1 ELECTRIC  FIELD  MEASUREMENTS 


Task  01 

Martin  Misakian 
Electricity  Division 

National  Institute  of  Standards  and  Technology 


1.1  Introduction 

The  objectives  of  this  project  are  to  develop  methods  to  evaluate  and  calibrate 
instruments  which  are  used,  or  are  being  developed,  to  characterize  the  electrical 
parameters  in  the  vicinity  of  power  lines  and  in  laboratory  apparatus  designed  to 
simulate  the  power  line  environment.  Electrical  measurement  support  is  also  provided 
for  the  Department  of  Energy  funded  elTorts  to  determine  if  there  are  biological  effects 
due  to  ac  fields  and  dc  fields  with  ions. 

During  1990,  NIST  organized  and  cohosted  with  Battelle  Pacific  Northwest 
Laboratories  (BPNL)  a workshop  on  exposure  and  biological  parameters  that  should 
be  considered  during  in  vitro  studies  with  extremely  low  frequency  (ELF)  magnetic 
and  electric  fields.  The  workshop  was  an  outgrowth  of  discussions  with  BPNL 
following  an  NIST  conference  presentation  late  in  1989.  It  was  evident  at  that 
time  that  all  of  the  relevant  exposure  parameters  during  many  in  vitro  bioeffect 
studies  with  ELF  magnetic  fields  were  not  being  adequately  characterized.  The 
NIST  contribution  to  the  workshop  considered  the  electric  field  and  current  density 
distributions  induced  in  culture  medium.  A portion  of  the  workshop  presentation  is 
reproduced  below.  Also  a related  NIST  paper  describing  an  optimum  experimental 
configuration  for  in  vitro  studies  with  ELF  magnetic  fields  was  published  [1]. 

In  addition  efforts  were  also  made  to  develop  equations  that  would  predict  the 
magnetic  field  in  a parallel  plate  magnetic  field  exposure  system  similar  in  principle 
to  one  developed  at  BPNL  [2].  The  results  of  the  calculations  are  described  below. 

Other  events  and  activities  related  to  the  NIST  Fields  Project  during  1990  included 
(a)  the  approval  and  publication  of  an  IEEE  standard  which  provides  guidance  for 
measuring  dc  electric  fields  and  ion  related  parameters  (IEEE  Std  1227-1990,  IEEE 
Guide  for  the  Measurement  of  DC  Electric-Field  Strength  and  Ion  Related  Quantities; 
draft  standard  prepared  at  NIST  in  collaboration  with  members  of  IEEE  DC 
Fields  and  Ions  Working  Group),  and  (b)  consultations  with  the  Food  and  Drug 
Administration  (electric  and  magnetic  field  measurements  near  electric  blankets), 
the  National  Toxicology  Program  (experimental  design  and  exposure  parameters 
for  rodent  study  with  magnetic  fields),  and  the  Environmental  Protection  Agency 
(reviews  of  draft  reports).  NIST  also  coauthored  two  working  group  papers  in  the 
IEEE  Power  Engineering  Society  related  to  characterizing  power  line  magnetic  fields. 
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Both  papers  are  scheduled  for  publication  in  the  IEEE  Transactions  on  Power  Delivery 
in  early  1991.  Site  visits  were  made  to  Southwest  Research  Institute  on  behalf  of  DOE 
and  to  the  University  of  Rochester  and  BPNL  on  behalf  of  EPRI,  during  which  electric 
and  magnetic  fields  were  characterized  in  biological  exposure  systems  and  calibration 
apparatus.  Preliminary  measurements  of  power  frequency  magnetic  fields  in  and  near 
the  NIST- Boulder  daycare  facility  were  performed.  The  measurements  were  coupled 
with  model  calculations  to  better  characterize  the  sources  of  magnetic  fields. 

NIST  serves  on  advisory  panels  for  the  California  Department  of  Health  Services  and 
the  National  Cancer  Institute  (NCI),  and  review  panels  for  the  National  Institutes 
of  Health  (NIH).  During  1990,  NIST  participated  in  meetings  of  the  NCI  and  NIH 
panels. 


1.2  Examination  of  Experimental  Designs  for  In  Vitro 
Studies  Using  ELF  Magnetic  Fields 

There  have  been  many  laboratory  studies  since  the  mid  1970’s  examining  the 
possibility  of  biological  effects  to  cells  in  vitro  due  to  power  frequency  magnetic 
fields.  These  studies  have  typically  used  Helmholtz  coils  to  produce  an  approximately 
uniform  magnetic  field  for  exposure  purposes.  Significantly,  the  enclosures  containing 
the  cells  and  liquid  growth  medium  have  varied  in  shape  and  size.  Because  the 
geometry  of  the  liquid  volume,  which  is  determined  by  the  enclosure,  has  a major 
impact  on  the  exposure  parameters,  the  enclosure  must  be  considered  as  an  important 
part  of  the  exposure  system.  This  portion  of  the  report  examines  the  exposure 
parameters  for  several  experimental  configurations  that  have  been  used  during  in  vitro 
studies  and  notes  some  of  their  advantages  and  limitations.  The  numerous  biological 
constraints  that  must  be  considered  for  the  different  configurations  will  be  the  subject 
of  a future  report.  Extensive  use  is  made  of  the  equations  developed  by  McLeod  et  al. 
[3]  for  determining  the  current  density  and  electric  field  induced  in  the  liquid  growth 
medium  by  the  magnetic  field. 

1.3  Experimental  Geometry  and  Exposure  Parameters 

In  vitro  studies  with  power  frequency  magnetic  fields  have  often  been  conducted  using 
circular  plastic  or  glass  petri  dishes  and  rectangularly  shaped  vessels  as  the  enclosures. 
When  the  enclosure  containing  the  liquid  and  biological  cells  is  introduced  into  the 
magnetic  field,  an  electric  field  and  an  associated  current  are  induced  in  the  liquid 
by  the  time-varying  magnetic  field.  Thus,  the  cells  are  exposed  to  three  candidate 
parameters  that  may  individually,  or  in  combination,  cause  a biological  effect,  i.e., 
the  magnetic  field  or  fiux  density,  B,  the  induced  electric  field  strength,  E,  and  the 
induced  current  density,  J . While  all  the  cells  will  experience  essentially  the  same 
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Figure  1.  Experimental  configurations  for  in  vitro  studies  with  ELF  magnetic  fields. 


magnetic  field,  the  exposure  to  E and  J will  depend  on  the  size  and  shape  of  the 
liquid  volume,  the  direction  of  the  magnetic  field  with  respect  to  the  licpiid  volume, 
and  the  location  of  the  cell  in  the  liquid  volume. 

We  consider  the  three  experimental  configurations  shown  in  figure  la:  a cylindrical 
liquid  volume  with  the  power  frequency  magnetic  field  parallel  and  perpendicular  to 
the  cylinder  axis  (figure  la  and  Ic)  and  a rectangular  liquid  volume  with  the  magnetic 
field  perpendicular  to  one  of  the  faces  (figure  lb).  It  should  be  noted  that  the  direction 
and  magnitude  of  any  dc  magnetic  field  may  be  of  significance  if  tests  of  a resonance 
model  are  being  conducted.  However  such  experiments  are  not  considered  here.  As 
qualitatively  indicated  in  figure  1,  the  induced  electric  field  lines  and  current  paths 
are  circular  when  the  magnetic  field  is  parallel  to  the  cylinder  axis.  In  general,  the 
electric  field  lines  and  current  paths  form  noncircular  loops  for  the  remaining  two 
cases.  The  equations  for  the  induced  current  density  and  electric  field  are  simplest 
when  the  magnetic  field  is  parallel  to  the  cylinder  axis.  Then  the  current  density  is 
given  by  [3] 
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J = (777  fBr,  (1) 

where  <7  is  the  conductivity,  / is  the  frequency  and  r is  the  radial  distance  from  the 
cylinder  axis.  Because  the  current  density  is  equal  to  the  product  crE^  the  electric 
field  strength  is  just 


E = JI(7  = 77  fBr.  (2) 

The  expressions  for  the  induced  current  density  in  the  rectangular  volume  are  more 
complicated  and  are  given  by  [1] 


and 
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where  2a  is  the  width  of  the  liquid  volume,  2h  is  the  depth,  and  u;  is  277  times 
the  frequency,  /.  The  ± sign  indicates  that  the  terms  in  the  summation  are 
alternatively  added  and  subtracted.  As  above,  the  corresponding  expressions  for 
the  electric  field  are  found  by  dividing  equations  (3)  and  (4)  by  the  conductivity, 
(7.  As  shown  by  McLeod  et  al.  [3],  the  case  of  the  cylindrical  volume  with  the 
magnetic  field  perpendicular  to  the  cylinder  axis  can  be  treated  by  cutting  the 
cylinder  into  approximately  rectangular  volumes  as  shown  in  figure  Ic  and  using 
equations  (3)  and  (4)  for  each  rectangular  volume.  The  rectangular  geometry  is 
examined  first. 


1.3.1  Rectangular  Geometry 

With  the  rectangular  configuration  shown  in  figure  lb,  it  can  be  shown,  using 
equations  (3)  and  (4),  that  the  induced  current  and  electric  field  vanish  along  the 
x-axis  as  well  as  along  the  four  edges  parallel  to  the  x-axis.  The  maximum  values 
of  the  induced  electric  field  and  current  density  occur  at  the  liquid  boundary  but 
the  locations  depend  on  the  cross  sectional  shape  of  the  volume.  The  cells  will  all 
experience  the  same  magnetic  field  because  the  magnetic  field  is  uniform  throughout 
the  liquid  volume.  However,  the  electric  field  and  current  density  experienced  by  the 
cells  will  vary  from  zero  to  some  maximum  value  depending  on  the  location  of  the 
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Figure  2.  Normalized  values  and  directions  of  F or  J along  bottom  surface  of  culture 
medium  when  B is  parallel  to  x-axis  (a)  culture  medium  is  2 mm  deep  (b)  culture  medium 
is  3 mm  deep. 

cell.  If  the  approximate  region  the  cells  occupied  during  exposure  to  the  magnetic 
field  were  known,  their  exposure  to  a more  limited  range  of  E and  J values  could  be 
determined  using  equations  (3)  and  (4).  Such  information  is,  in  general,  not  known 
when  experiments  are  conducted  with  cells  in  suspension. 

Our  ability  to  describe  the  exposure  conditions  for  the  cells  improves  if  the  cells  can 
be  cultured  on  the  bottom  surface  of  the  enclosure.  Then  the  vertical  components 
of  E and  J vanish  and  approximately  uniform  values  of  E and  J in  the  horizontal 
direction  can  be  produced  along  the  bottom  surface  of  the  liquid  if  the  liquid  depth 
is  small  compared  to  the  width,  i.e.,  2h  2a.  Figure  2a  shows  normalized  values 
of  FJ  or  J calculated  at  the  bottom  surface  of  the  liquid  volume  for  2a  = 6 cm 
and  2h  = 0.2  cm.  It  can  be  seen  for  this  example  that  the  region  of  approximately 
uniform  E oi  J extends  almost  to  the  side  edges  where  E and  J rapidly  decrease 
to  zero.  Increasing  the  depth  of  the  liquid  by  1 mm  decreases  the  uniformity,  but 
because  the  liquid  depth,  2/i,  is  still  much  smaller  than  the  side  dimension,  2a,  the 
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uniformity  is  still  very  good  as  shown  in  figure  2b.  While  the  1 mm  difference  in  depth 
has  had  a small  effect  on  the  uniformity  of  E and  J in  this  example,  the  increase  in 
magnitude  of  E and  J is  not  insignificant,  amounting  to  50  % over  most  of  the  area. 
This  sensitivity  of  E and  J to  liquid  depth  for  the  rectangular  configuration  may  be 
an  important  consideration  when  attempts  are  made  to  replicate  the  results  of  an 
earher  in  vitro  study. 

It  should  also  be  noted  that  cells  cultured  on  the  bottom  surface  of  the  enclosure 
will  experience  an  “uneven”  exposure  to  the  approximately  uniform  E and  J values 
described  above  because  a portion  of  the  underside  of  the  cell,  which  is  in  contact 
with  plastic  or  glass,  will  not  be  exposed  to  the  electric  field  or  current  density  in  the 
liquid.  This  problem  can  be  reduced  if  the  cells  can  be  cultured  in  agar  that  has  the 
same  conductivity  as  the  liquid.  However,  the  use  of  agar  can  introduce  undesirable 
biological  effects  [4]  and  cause  significant  changes  in  the  exposure  conditions.  For 
example,  if  liquid  2 mm  deep  covers  a layer  of  agar  2 mm  thick  on  the  bottom  of  a 
rectangular  enclosure  3 cm  wide,  and  if  cells  are  cultured  near  the  upper  surface  of 
the  agar,  they  will  experience  little  or  no  exposure  to  £*  or  J except  near  the  side 
edges. 

Summarizing  _ the  above  discussion,  the  rectangular  configuration  allows  for  the 
exposure  of  large  numbers  of  cells  to  uniform  values  of  B,  B,  and  J if  they  are 
cultured  on  the  bottom  surface  of  the  enclosure.  The  values  of  E and  J are  sensitive 
to  liquid  depth  and  can  be  calculated  using  equations  (3)  and  (4).  The  equations 
for  determining  E and  J are  not  elementary,  but  can  be  readily  evaluated  with  a 
computer.  The  use  of  agar  can  significantly  influence  the  exposure  parameters  can 
and  also  have  biological  consequences. 

1.3.2  Cylindrical  Geometry:  Magnetic  Field  Perpendicular  to  Axis 

As  noted  earlier,  the  case  of  a cylindrical  liquid  volume  with  the  magnetic  field 
perpendicular  to  the  cylindrical  axis  can  be  treated  by  imagining  the  cylinder  being 
cut  into  rectangular  sections  and  using  the  theory  for  rectangular  geometry.  Thus, 
the  above  remarks  regarding  the  characterization  of  exposure  conditions  for  cells  in 
suspension  and  cells  cultured  on  the  bottom  of  the  enclosure  again  apply  to  each 
section.  Figure  3 shows  normalized  values  of  E and  J at  the  bottom  surface  of  a 
liquid  volume  6 cm  in  diameter  and  0.2  cm  deep.  The  uniformity  is  again  seen  to  be 
very  good  over  most  of  the  surface  area  and  is  a consequence  of  the  fact  that  2h  is 
still  much  smaller  than  the  width  of  the  imaginary  sections  at  most  locations.  As  for 
the  examples  considered  in  figure  2,  the  magnitude  of  B or  J is  sensitive  to  the  liquid 
depth. 


1.3  Experimental  Geometry  and  Exposure  Parameters 
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Normalized  E,  J Values  at  Bottom 
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Figure  3.  Normalized  values  and  directions  of  E and  J along  bottom  surface  of  culture 
medium  when  B is  parallel  to  x-ajcis. 
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Normalized  E,  J Values  in 
Cylindrical  Volume  of  Liquid 

diameter  = 6 cm 
B II  to  cylinder  axis 

Figure  4.  Normalized  values  and  directions  of  E and  J in  culture  medium  when  B is 
parallel  to  z-axis 

1.3.3  Cylindrical  Geometry:  Magnetic  Field  Parallel  to  Axis 

When  the  magnetic  field  is  parallel  to  a cylindrical  volume,  the  problem  of  describing 
the  exposure  conditions  for  the  cells  becomes  more  simple.  The  equations  describing 
E and  J (equations  (1)  and  (2))  are  elementary  and  indicate  that  both  E and  J 
increase  linearly  with  radial  distance  from  the  cylinder  axis.  As  noted  earlier,  the 
electric  field  lines  and  current  density  form  concentric  circular  loops  and  are  thus  more 
easily  visualized  than  for  the  rectangular  geometry.  Significantly,  equations  (1)  and  (2) 
indicate  that  E and  J are  independent  of  liquid  or  agar  depth. 

Figure  4 shows  normalized  E and  J values  for  a cylindrical  volume  6 cm  in  diameter. 
While  the  lack  of  uniformity  of  these  exposure  parameters  may  at  first  seem  to  be 
a disadvantage  when  characterizing  cell  exposure,  it  is  possible  to  make  use  of  the 
nonuniformity  to  obtain  several  ranges  of  exposure  to  E and  J during  the  same 
experimental  run  [1].  For  example,  cells  can  be  segregated  into  annular  sections  with 
circular  walls  and  the  range  of  E and  J values  experienced  by  the  cells  in  each  section 
can  be  determined  from  equations  (1)  and  (2)  and  the  radii  that  define  the  section. 
The  exposure  geometry  and  design  of  an  enclosure  that  would  be  suitable  for  short 
term  studies  with  cells  in  suspension  are  shown  in  figure  5. 

While  the  cylindrical  volume  in  figure  5 has  been  partitioned  into  three  sections,  more 
or  fewer  sections  could  be  used.  For  more  extended  periods  of  exposure,  lasting  as 


1.3  Experimental  Geometry  and  Exposure  Parameters 
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Figure  5.  Experimental  design  for  short  term  in  vitro  studies  with  cells  in  suspension. 
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Annular  regions  indicated  with  circles 
painted  on  bottom  of  petri  dish 


Figure  6.  Method  for  characterizing  radial  positions  of  cells. 

long  as  six  or  seven  days,  the  enclosure  can  be  fitted  with  a cap  and  “sealed”  with 
commercially  available  laboratory  film  which  prevents  significant  evaporation  but 
allows  adequate  air  exchange  [5].  Vertical  stacking  of  the  enclosures  is  one  technique 
for  increasing  the  number  of  cells  that  are  exposed  during  a single  experimental  run. 
Introducing  a layer  of  agar  at  the  bottom  of  the  enclosure  as  shown  in  figure  5 provides 
for  a more  uniform  exposure  of  the  cells  (top  and  underside)  to  E and  J . Unlike 
the  rectangular  geometry  case,  introduction  of  the  agar  will  not  affect  the  values  of 
E and  J because  E and  J are  independent  of  liquid  or  agar  depth.  However,  use  of 
agar  can  have  significant  biological  effects  as  noted  previously. 

If  the  circular  walls  shown  in  figure  5 introduce  an  impediment  for  cells  plated 
directly  on  the  enclosure  surface,  it  may  be  possible  to  replace  the  walls  with  painted 
concentric  circles  on  the  underside  of  the  enclosure  as  shown  in  figure  6. 

Other  techniques  for  characterizing  the  radial  positions  of  the  ceils  in  the  cylindrical 
enclosure  include  growing  the  cells  on  a glass  cover  slip  placed  in  the  enclosure 
(the  cover  slip  can  be  removed  later  for  measurements  with  a rule  and  compass) 
or  indicating  the  annular  regions  on  a flat  surface,  as  shown  in  figure  7,  upon  which 
the  enclosure  is  placed. 

Because  of  the  spatial  uniformity  of  B and  the  readily  determined  variation  of  E and  J 
with  position  in  the  enclosure,  the  experimental  configuration  described  in  this  section 
can  be  used  to  investigate  directly  questions  regarding  the  relative  significance  of 
interaction  mechanisms  involving  the  magnetic  field  and/or  induced  electric  field  and 
current  density.  For  example,  the  observation  of  a biological  effect  at  the  same  level 
in  all  sections  of  the  enclosure  would  suggest  a mechanism  involving  direct  magnetic 
interaction,  whereas  an  effect  that  differed  between  sections  could  be  indicative  of  a 
mechanism  involving  the  induced  electric  field  and  magnetic  field  combination.  For 
the  latter  case,  a single  experiment  could  provide  data  for  a range  oi  E oi  J exposure 
levels. 


1.4  Magnetic  Fields  from  Parallel  Plates 
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Annular  regions  indicated  with 
circles  painted  in  well 


Figure  7.  Method  for  characterizing  radial  positions  of  cells. 

Summarizing  this  section,  the  cylindrical  configuration  with  the  magnetic  field  parallel 
to  the  cylinder  axis  results  in  induced  electric  fields  and  current  densities  that 
are  easily  visualized,  and  elementary  equations  for  calculating  E and  J that  are 
independent  of  liquid  depth.  In  addition,  the  amount  of  information  pertaining 
to  dose-response  during  a single  experimental  run  can  be  maximized,  and  it  may 
be  possible  to  distinguish  between  a purely  magnetic  field  bioeffect  and  a bioeffect 
that  involves  the  electric  field  or  current  density-magnetic  field  combination.  The 
configuration  in  this  section  does  not  permit  exposure  of  large  numbers  of  cells  to  the 
same  E and  J values.  The  use  of  agar  does  not  have  a significant  effect  on  E and  J 
values,  but  can  have  biological  effects  as  noted  earlier. 

1.4  Magnetic  Fields  from  Parallel  Plates 

Magnetic  fields  for  ELF  bioeffect  studies  have  typically  been  produced  with  solenoids, 
Helmholtz  coils,  or  rectangular  coils  with  a common  axis.  Because  stray  fields  from 
such  coil  systems  can  extend  to  considerable  distances  and  because  of  the  difficulty  in 
shielding  sham  exposed  animals  or  cells  from  these  fields,  it  is  of  interest  to  consider 
an  experimental  apparatus  recently  reported  by  Miller  et  al.  [2],  that  has  low  stray 
fields.  Figure  8 shows  a schematic  view  of  the  apparatus  which  consists,  in  part, 
of  three  stacked  parallel  plates.  The  ac  current  in  the  central  plate  divides  equally 
between  the  upper  and  lower  plates  and  produces  a magnetic  field  parallel  to  the 
plates  in  each  half  of  the  apparatus  as  shown  schematically  in  figure  8 at  some  instant. 
Because  the  fields  in  the  two  sections  are  in  opposite  directions,  stray  fields  from  the 
apparatus  decrease  more  rapidly  than  the  field  from  either  section  alone.  This  section 
describes  attempts  to  theoretically  describe  the  magnetic  field  for  a similar  parallel 
plate  apparatus. 
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Figure  8.  Schematic  view  of  parallel  plate  apparatus  (after  Ref.  2)  showing  direction  of  ac 
current  at  some  instant  and  direction  of  magnetic  held  between  parallel  plates. 


Figure  9.  Half-section  view  of  current  distribution  in  hat  conductor  that  changes  in  width 
by  a factor  of  four. 

Efforts  to  calculate  the  held  for  the  apparatus  shown  in  hgure  8 are  complicated  by  the 
nonuniformity  of  the  current  density  at  the  points  of  connection  between  the  plates 
and  the  electrodes  that  provide  current.  Figure  9 shows  how  the  current  distribution 
departs  from  uniformity  when  a single  connecting  electrode  is  one-fourth  as  wide  as 
the  plate  [6]. 

The  nonuniformity  of  the  current  density  due  to  connections  of  metal  of  uneven  widths 
can  be  overcome  by  the  modihcation  shown  in  hgure  10.  The  electrodes  providing 
current  to  the  plates  are  expanded  in  width  to  match  the  parallel  plates  and  the 
shape  has  been  changed  to  resemble  two  rectangular  solenoids.  Attempts  were  made 
to  determine  the  magnetic  held  theoretically  by  applying  the  law  of  Biot-Savart  to  just 
the  upper  half  of  the  arrangement  as  shown  in  hgure  11.  Because  of  the  quasistatic 
nature  of  the  held,  it  was  assumed  that  the  time  dependence  of  the  magnetic  held 


Figure  10.  Modihed  parallel  plate  apparatus. 


1.4  Ma,gnetic  Fields  from  Parallel  Plates 
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Figure  11.  Geometry  and  coordinates  for  upper  half  of  parallel  plate  exposure  apparatus. 
The  law  of  Biot-Savart  takes  the  form  of  a double  integral  (example  shown  above)  and  is 
applied  to  aU  surfaces  of  the  apparatus. 

could  be  predicted  by  multiplying  the  dc  magnetic  field  (obtained  from  the  law  of 
Biot-Savart)  by  a sincjt  term. 

Figure  12  shows  comparisons  of  measured  horizontal  magnetic  field  profiles  at 
frequencies  of  25  Hz  and  1 kHz  in  a rectangular  apparatus  as  shown  in  figure  11 
with  dimensions  of  50  cm  x 50  cm  x 10  cm  side  dimension.  The  measurements  were 
performed  with  a miniature  coil  probe  5.7  mm  long  and  4.4  mm  in  diameter. 
Agreement  is  fair  at  25  Hz  and  poor  at  1 kHz.  At  60  Hz,  the  magnetic  field  is  about 
3.5  % less  than  the  dc  value  at  the  center  of  the  apparatus.  The  departure  from  the 
dc  theory  is  due  to  a skin  effect  [7]  which  causes  the  current  to  concentrate  along  the 
edges  of  the  plate  as  shown  schematically  in  figure  13. 

The  concentration  results  because  of  the  larger  inductance  (and  correspondingly 
increased  impedance)  along  a central  path  through  the  plates.  It  is  noted  that  even 
with  nonuniform  current  distributions,  regions  of  approximately  uniform  magnetic 
field  can  be  obtained  between  the  parallel  plates  with  an  apparatus  similar  to  that 
shown  in  figure  8 [2].  The  field  magnitude  and  uniformity  must  be  determined  by 
measurement,  however. 

The  problems  with  current  nonuniformity  (including  the  skin  effect)  can  be  overcome 
by  replacing  the  metal  plates  with  two  adjacent  rectangular  solenoids  as  shown  in 
figure  14. 
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Figure  12.  Horizontal  magnetic  field  profiles  midway  between  top  and  bottom  plates  of 
parallel  plate  apparatus  50  cm  x 50  cm  x 10  cm  spacing.  Magnetic  field  values  for  dc,  25  Hz, 
and  1 kHz  frequencies  are  indicated. 


Figure  13.  Edge-on  view  of  plate  showing  increased  current  density  concentration  at  edges 
of  plate. 


Figure  14.  Arrangement  of  “coupled”  rectangular  solenoids  to  reduce  stray  magnetic  fields. 


i.4  Magnetic  Fields  from  Parallel  Plates 


The  fields  in  each  section  as  well  as  the  stray  fields  from  the  '‘coupled”  arrangement 
can  be  calculated  using  equations  for  rectangular  loops  of  wire  [8]  and  the  principle 
of  superposition.  The  design  shown  in  figure  14  eliminates  significant  skin  effects 
and  requires  less  current  than  the  plates  to  produce  the  same  field  because  of  the 
many  turns  of  wire.  As  for  the  parallel  plate  arrangement,  the  stray  field  from  the 
“coupled”  rectangular  solenoids  will  decrease  faster  than  either  solenoid  alone. 
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2.1  Introduction 

The  objectives  of  this  project  are  the  development  of  measurement  capabilities  and 
the  providing  of  fundamental  data  as  part  of  the  Department  of  Energy’s  basic 
research  concerned  with  the  development  and  evaluation  of  advanced  compressed- 
gas-insulation  technology. 

To  reduce  space  requirements  and  improve  the  efficiency  of  high-voltage  transmission 
systems,  the  electric  power  industry  has  turned  toward  more  extensive  use  of 
compressed-gas  insulation.  To  design  meaningful  tests  of  system  performance  and 
establish  specifications  for  the  quality  of  materials  used  in  such  systems,  more 
information  is  needed  about  the  fundamental  physical  and  chemical  processes 
which  lead  to  insulation  deterioration  and  electrical  breakdown.  The  research 
includes  applications  of  gas  chromatography-mass  spectrometry  to  characterize 
corona  discharge  by-products;  and  the  acquisition  of  fundamental  data,  such  as, 
reaction-rate  coefficients,  corona-inception  voltages,  production  rates  of  corona  by- 
products, the  effects  of  contaminants  on  discharge  initiation,  and  the  rates  of 
discharge-induced  decomposition  of  the  gas. 

This  section  of  the  report  highlights  the  progress  of  research  in  two  areas.  The  first 
is  concerned  with  the  development  of  a sensitive  method  for  detecting  trace  levels 
of  disulfur  decafluoride  (S2F10)  in  sulfur  hexafluoride  (SFe)  by  gas  chromatography 
mass  spectrometry.  The  second  area  of  research  has  to  do  with  development  of  an 
improved  stochastic  analyzer  for  pulsating  phenomena  (SAPP)  and  its  application  to 
the  investigation  of  negative-corona  discharge  pulses  in  gases. 
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The  work  on  S2F10  detection  was  carried  out  in  collaboration  with  scientists  at 
the  Oak  Ridge  National  Laboratory  and  results  have  already  been  reported  in  one 
archival  [9]  and  two  conference  papers  [10, 11].  The  work  on  stochastic  behavior  of 
corona  pulses  involved  participation  by  guest  scientists  from  the  Plasma  Research 
Institute  of  India  and  Old  Dominion  University.  Results  from  the  latter  investigation 
have  appeared  in  one  archival  [12]  and  two  conference  papers  [13, 14].  Three  other 
archival  papers  covering  different  aspects  of  this  work  have  also  recently  been 
submitted  for  publication  [15-17]. 


2.2  Detection  of  Trace  S2F10  in  SFe 

The  ability  to  detect  gaseous  disulfur  decafluoride  (S2F10)  at  trace  levels  in  sulfur 
hexafluoride  (SFs)  has  become  of  interest  because  of  its  known  high  level  of 
toxicity  [18-22]  and  because  of  recent  measurements  [23,  24]  and  calculations  [25,  26] 
which  show  that  S2F10  can  be  produced  during  electrical-discharge  activity  in  SFe. 
The  concern  about  S2F10  production  in  SFs  stems  in  part  from  the  increasingly  wide 
spread  use  of  SFe  as  a gaseous  dielectric  in  high-voltage  electric-power  transmission 
'systems.  The  Occupational  Safety  and  Health  Administration  (OSHA)  has  recently 
set  peak  exposure  limits  (PEL)  for  S2F10  at  lOppb  [27].  However,  due  to  the  difficulty 
of  detecting  low  concentrations  of  S2F10  in  SFe,  enforcement  of  the  PEL  has  been 
stayed  until  the  development  of  suitable  detection  techniques. 

There  are  previous  preliminary  reports  by  Janssen  [28, 29]  of  having  achieved  a 
detection  sensitivity  for  S2F10  in  helium  gas  of  0.04  ppb  by  using  a gas  chromatograph 
equipped  with  an  electron  capture  detector.  In  order  to  achieve  this  level  of  sensitivity, 
however,  it  was  necessary  to  employ  a complex  S2F10  enrichment  procedure  requiring 
long  analysis  times  (~  45  minutes).  When  SFe  is  a major  constituent  of  the  gas  in 
which  S2F10  is  to  be  detected,  the  sensitivity  for  S2F10  detection  drops  dramatically 
(to  > 10  ppb)  despite  the  enrichment  process.  This  drop  in  sensitivity  is  a consequence 
of  increased  background  signal  in  the  detector  output  due  to  the  favorable  electron 
attaching  properties  of  the  residual  SFe  which  interferes  with  the  S2F10  signal. 
Reduced  sensitivity  for  S2F10  detection  in  the  presence  of  SFe  also  occurs  when  flame 
photometric  detectors  [28]  or  thermal  conductivity  detectors  [11]  are  used. 

Mass-spectrometric  detection,  when  coupled  to  a gas  chromatograph,  suffers  from  the 
specific  limitation  that  S2F10  has  a positive-ion  mass  spectrum  which  is  very  similar 
to  that  of  SFe,  especially  at  the  standard  electron-impact  energy  of  70  eV  [29,  30,  31]. 
For  electron-impact  energies  greater  than  20  eV,  all  of  the  significant  ions  that  appear 
in  the  S2F10  mass  spectrum  also  appear  in  the  SFe  spectrum.  Thus,  even  when 
chromatographic  separation  of  S2F10  from  SFe  is  used,  ions  from  SFe  can  sufficiently 
increase  the  background  level  to  interfere  with  S2F10  detection.  Recent  claims  [32]  of 
using  a mass  spectrometer  operated  at  low  electron-impact  energies  to  detect  trace 
levels  of  S2F10  in  SFe  have  been  disputed  [31]. 
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Figure  15.  Schematic  diagram  of  GC-MS  system  with  a heated  jet  separator  between  the 
GC  column  and  the  ion  source  of  the  mass  spectrometer.  Heating  of  the  interface  between 
the  GC  and  MS  is  essentiiil  for  the  sensitive  detection  of  S2Fi()  in  the  presence  of  SF^. 

.A.  new  method  is  described  here  in  which  a gas  chromatograph-mass  spectrometer 
(GC-MS)  system  is  easily  modified  to  allow  detection  of  S2Fio  in  gaseous  SF^  down  to 
the  ppb  level  without  the  necessity  of  a complex  and  time-consuming  pre-enrichment 
procedure.  With  this  method,  a fraction  of  the  S2Fio  is  converted  to  thionyl  fluoride 
(SOF2)  in  a heated  jet-separator  at  the  interface  between  the  GC  and  MS.  Because  the 
major  ions  that  appear  in  the  70 -eV  mass  spectrum  of  SOF2  differ  from  those  in  the 
SFfi  spectrum,  interference  of  SF^  with  the  S2F10  peaks  in  single-ion  chromatograms 
can  be  avoided.  Elimination  of  the  enrichment  step  significantly  reduces  the  gas 
analysis  time  and  improves  the  reliability  of  the  detection  method. 

2.2.1  Measurern,tnt  System 

The  analytical  method  proposed  for  sensitive  detection  of  S2F10  in  SFg  utilizes  a 
GC-MS  with  a heated  metal  tube  on  the  lower-pressure  side  of  the  interface  between 
the  GC  and  MS.  As  shown  in  figure  15,  a jet  separator  serves  as  the  heated  metal  tube 
in  the  GC-MS  system  used  to  obtain  the  residts  presented  here.  The  jet  separator  is  a 
standard  design  consisting  of  a stainless  steel  tube  of  8 cm  length  and  5 mm  diameter 
which  is  heated  to  a temperature  of  about  180 °C.  As  discussed  later,  a fraction 
of  the  S2F10  undergoes  thermal  decomposition  in  the  heated  tube  and  converts  to 
SOF2  via  surface-catalyzed  wall  reactions  with  H2O.  The  SOF2  can  then  be  detected 
by  the  mass  spectrometer  without  interference  by  the  ion  signals  from  SFg.  The 
efficiency  of  the  conversion  is  not  known  but  appears,  from  experiments  discussed 
below,  to  depend  on  the  length,  shape,  and  temperature  of  the  heated  surface  upon 
which  conversion  takes  place.  The  conversion  efficiency  is  most  sensitive  to  the  tube 
temperature  and  relatively  insensitive  to  details  of  the  tube  geometry  and  material 
composition.  A high  sensitivity  for  S2F10  detection  is  obtainable  with  a wide  range 
of  heated  inlet  tube  configurations. 
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The  GC-MS  system  used  for  the  experiments  presented  here  is  a Hewlett  Packaxd’ 
5992A  with  a Hewlett  Packard  jet  separation  inlet.  In  this  particular  instrument 
the  jet  separator  is  heated  by  the  diffusion  pump  heater.  For  other  systems  the 
conversional  region  may  be  heated  externally.  The  GC  column  is  similar  to  that 
described  by  Hanrahan  and  Patterson  [33],  namely  a 24'  x 1/8"  Teflon  tube  containing 
30  % SP-2100  fSupelcoj  on  80/100  chromosorb  WAW.  The  column  was  operated  using 
helium  as  a carrier  gas  (20-30  ml/min.  flow  ratej  which  had  not  been  predried.  The 
GC  oven  temperature  was  maintained  between  25  and  50  '^C  and  the  sample  injections 
were  made  using  a gas-tight  syringe.  Detection  sensitivities  for  S2F]o  in  SF^  of  nearly 
1 ppb  were  obtainable  using  the  instrument  described  above. 

Pure  S2F10  samples,  synthesized  at  Clemson  University,  were  provided  to  us  as  liquid 
under  its  own  vapor  pressure  (~  600  torr  at  room  temperature j in  30  ml  stainless  steel 
cylinders.  Identity  and  purity  of  the  sample  were  established  by  IR  spectroscopy  and 
GC-MS  analysis.  Gas-phase  reference  mixtures  were  prepared  immediately  prior  to 
use,  either  by  injection  or  needle  valve,  in  argon  (99.999  %j  or  SFq  (99.99  %J  in  a 
stainless  steel  vacuum  manifold  housed  in  a vented  hood.  The  argon  and  SF^  were 
not  dried  before  use. 


2.2.2  Data  Acquisition  and  Analy.si3 

Data  acquisition  was  performed  using  the  manufacturer-supplied  software  and 
analysis  was  performed  using  software  developed  in-house.  For  S2F10  samples  in 
buffer  gcLses  other  than  SFq,  trace  amounts  of  S2F]o  could  be  detected  by  single  ion 
monitoring  of  SFf , SF^,  SFf , SF^  and  F'^  (rn/ z — 127,  108,  89,  70,  and  19j  produced 
by  direct  ionization  of  S2Fio[31].  When  mixed  with  SF^,  sensitive  S2F]o  detection 
required  the  monitoring  of  ion  signals  characteristic  of  the  SOF2  produced  by  surface 
conversion  in  the  heated  inlet  tube  mentioned  above.  These  ions  include  SOF-]", 
SOF"^,  or  (m/z  = 86,  67,  and  48j.  For  the  data  presented  here,  as  many  as  8 
ions  could  be  monitored  at  a given  time  with  dwell  times  ranging  from  10  to  200  ms. 

The  response  of  the  GC-MS  to  S2F10  is  defined  here  to  be  the  area  under  the  peak 
identified  with  this  species  in  a chromatogram.  This  area  was  computed  from  the 
GC-MS  data  by  numerical  integration  after  performing  a background  subtraction. 
No  attempts  were  made  to  “smooth”  or  average  data  and  unless  specified,  the 
data  presented  here  corresponds  to  chromatograms  from  single  sample  injections. 
Retention  time  was  measured  from  the  time  of  sample  injection. 

^The  identification  of  commercial  materials  and  their  sources  is  made  to  describe  the  experiment 
adequately.  In  no  ca.se  does  this  identification  imply  recommendation  by  the  National  Institute  of 
Standards  and  Technology,  nor  does  it  imply  that  the  product  is  the  best  available. 
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Figure  16.  Chromatograms  of  representative  ion  signals  due  to  ions  from  direct  ionization 
of  S2F  10  s-ud  from  ionization  of  surface  conversion  products  for  a 500  ppm  sample  of  S2F10 
in  100  kPa  argon  buffer  gas. 

2.2.3  GC-MS  Results  and  Calibration 

When  a gas  sample  containing  S2F10  is  injected  into  the  column  of  the  GC-MS  in 
which  the  interface  is  heated,  the  conversion  of  S2F10  into  SOF2  discussed  in  the 
next  section  becomes  evident  from  single-ion  chromatograms  at  mj z = 48,  67,  and  86 
(SO"^,  SOF'^,  and  SOFj).  Examples  of  such  chromatograms  are  shown  in  figure  16 
for  a 500  ppm  sample  of  S2F10  in  argon.  Ions  produced  by  direct  ionization  of  S2F10 
are  observed  (F'*’  and  SF^)  as  are  ions  from  SOF2  produced  by  surface  conversion 
(SOF+  and  SO+). 

The  mjz  = 20  chromatogram  in  figure  16  also  exhibits  a peak  at  the  retention 
time  characteristic  for  S2F10,  thus  indicating  that  HF  is  also  a by-product  of  S2F10 
decomposition  in  the  heated  tube.  This  peak,  however,  exhibits  structure  which  can 
be  accounted  for  by  interference  from  HF  originally  present  as  a contaminant  in  the 
sample  mixture.  Analysis  of  dilute  HF  mixtures  showed  that  HF  has  nearly  the  same 
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Figure  17.  Single  ion  cliromatograms  of  a 660  ppb  S2F10  sample  in  200  kPa  argon  buffer 
gas. 


retention  time  as  S2F10  for  the  column  conditions  used  here.  Thus,  although  HF  may 
be  formed  by  thermal  decomposition  of  S2F10  in  the  heated  gas  inlet  tube,  it  cannot, 
by  itself,  be  relied  upon  to  give  an  unambiguous  indication  of  S2F10. 

For  the  results  shown  in  figure  16,  ions  due  to  direct  ionization  of  S2F10  are  observable 
because  no  SFe  was  present  in  the  sample.  However,  in  the  presence  of  SFs,  S2F10 
signals  in  the  single-ion  chromatograms  at  m/z  = 19  and  127  are  obscured  by  the  SFe 
contribution  to  the  signal  at  these  masses.  This  is  demonstrated  by  data  presented 
in  figure  17  from  a gas  sample  containing  660  ppb  S2F10  in  SFe.  These  results  show 
the  effect  of  SFe  in  limiting  the  sensitivity  for  S2F10  detection  when  conventional 
mass-spectrometric  monitoring  is  used,  and  clearly  show  the  advantages  of  using  the 
S2F10  conversion  technique  discussed  here. 

Figure  18  shows  a chromatogram  for  m/z  = 86  for  an  "unknown*'  SFe  sample  from  a 
container  in  which  the  gas  had  been  subjected  to  a weak  electrical  corona  discharge. 
The  discharge  produced  trace  amounts  of  the  ox}*huorides  SOF2,  SOF4,  and  SO2F2  in 
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Figure  18.  Single  ion  chromatogram  of  an  “unknown”  sample  of  SFe  that  was  partially 
decomposed  in  a dc  corona  discharge. 
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Figure  19.  Detection  response  curves  for  various  concentrations  of  S2F10  and  SOF2  in 
argon  using  a chromosorb  GC  column  discussed  in  the  text.^  Open  symbols  represent 
detection  of  S2F10  and  SOF2  by  direct  ionization,  while  closed  symbols  represent  detection 
of  S2F10  by  ionization  of  SOF2  produced  by  surface  catalysis.  The  measured  responses  were 
normalized  to  the  observed  responses  of  a reference  sample  containing  750  ppm  S2F10  and 
750 ppm  SOF2  in  argon.  Similar  response  curves  are  obtained  for  samples  diluted  in  SFe. 

addition  to  S2F10.  It  is  seen  that  the  retention  times  of  the  SOF2  and  SOF4  produced 
by  the  discharge  are  much  shorter  than  for  S2F10  (2.1  and  2.3  minutes  compared  to 
4.0  minutes).  The  oxyhuorides  are  known  from  previous  work  [34-36]  to  be  the  major 
gaseous  by-products  from  decomposition  of  SFe  in  corona  discharges  when  at  least 
trace  levels  of  oxygen  and  water  vapor  are  present.  The  production  of  S2F10  in  SFe 
electrical  discharges  has  only  recently  been  verified  [23,24]. 

Under  the  conditions  used  to  obtain  the  data  in  figures  16  and  17,  the  response  of 
the  GC-MS  was  measured  as  a function  of  the  S2F10  concentration  in  reference  gas 
samples.  Typical  response  curves  obtained  using  known  concentrations  of  S2F10  mixed 
with  argon  are  shown  in  figure  19  for  the  indicated  m/^  values.  The  S2F10  response 
is  compared  with  that  for  SOF2  which  was  also  added  to  the  mixture  in  known 
concentrations.  The  S2F10  response  curves  for  miz  = 48  and  86,  corresponding  to 
the  chemical  conversion  of  S2F10  into  SOF2,  are  seen  to  be  distinctly  nonlinear.  The 
deviation  from  linearity  is  greatest  at  the  highest  S2F10  concentrations  and  indicates 
a saturation  effect  possibly  associated  with  depletion  of  H2O  from  the  heated  jet 
separator  tube  surface  during  the  conversion  of  S2F10  into  SOF2.  By  contrast,  the 
S2F10  response  curve  for  m/z  = 127,  corresponding  to  unreacted  S2F10,  is  seen  to  be 
nearly  linear  over  5 orders  of  magnitude.  Likewise  the  responses  to  SOF2  originally 
present  in  the  sample  at  m/z  = 48  and  36  are  also  nearly  linear.  The  small  deviations 
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from 'linearity  observed  at  low  concentrations  (<  200  ppb)  may  be  attributable  to 
uncertainties  in  the  preparation  of  samples  by  successive  dilutions. 

The  shape  and  slope  of  the  response  curves  for  the  S2F10  conversion  process  are  quite 
reproducible  for  a given  set  of  instrument  operating  conditions.  Nevertheless,  when 
attempting  accurate  measurements  of  S2F10  concentrations  in  SFe,  frequent  checks 
on  the  GC-MS  calibration  may  be  necessary  and  reference  sample  concentrations 
should  ideally  be  as  close  to  the  unknown  sample  concentration  as  possible.  The 
response  of  the  present  GC-MS  system  to  S2F10  was  periodically  checked  and  found 
to  be  stable  for  long  periods  of  time  (tiq^ically  more  than  a day).  WTien  performing  a 
quantitative  analysis  of  S2F10  content  in  SFe,  it  is  desirable  to  bracket  the  unknown 
sample  injection  with  reference  injections  of  comparable  concentration. 

It  should  be  noted  that  although  the  response  at  high  S2F10  concentrations  shows 
evidence  of  H2O  removal,  the  surface  apparently  recovers  to  its  former  condition 
in  a relatively  short  time  5-10  min.)  as  indicated  by  responses  from  successive 
injections  of  samples  at  high  concentrations  100  ppm).  In  general,  the  larger  the 
quantity  of  S2F10  injected,  the  longer  the  surface  recovery  time  of  the  inlet  tube.  For 
samples  containing  less  than  10  ppm  S2F10,  the  recovery  time  of  the  system  used  here 
was  found  to  be  less  than  the  typical  S2F10  retention  time  of  4 min. 

The  limit  of  S2F10  detection  in  SFe  depends,  of  course,  on  the  operating  conditions  of 
the  GC-MS  equipped  with  the  heated  inlet  line.  For  the  instrument  used  here,  under 
conditions  identical  to  those  used  to  obtain  the  results  shown  in  figures  16, 17  and  IS, 
the  limit  of  detection  (signal- to- noise  ratio  = 2)  was  estimated  to  be  about  2 ppb 
for  a 2 ml  injection  at  200  kPa.  Figure  20  shows  an  example  of  a result  obtained  for 
a 4 ppb  mixture  of  S2F10  in  SFs  which  gives  an  indication  of  the  limit  of  detection 
sensitivity  for  this  instrument.  The  single-ion  chromatogram  at  m/z  = S6  is  an 
average  of  data  from  four  injections  each  with  a volume  of  2 ml  at  a total  gas 
pressure  of  200  kPa.  The  S2F10  peak  is  statistically  significant  and  is  indicated  by 
the  vertical  arrow.  The  presence  of  S2F10  in  SFe  at  the  10 ppb  level  (the  PEL) 
gives  a much  larger  peak  than  that  shown  in  figure  20  and  is  easily  detected  using  a 
single  injection.  If  necessary,  additional  improvements  in  signal-to-noise  ratios  needed 
for  higher  sensitivity  applications  can  be  achieved  through  more  signal  processing, 
(i.e.,  averaging  of  results  from  more  injections  coupled  %vith  appropriate  background 
subtraction)  or  using  heated  MS  inlet-tube  conditions  that  optimize  the  efficiency  of 
S2F10  conversion  into  SOF2. 


2.2.4  M ass- Spectrometric  Investigation  of  S2F10  Pyrolysis 

In  order  to  determine  the  optimum  operating  temperature  range  for  the  heated 
gas  inlet  tube  and  to  better  understand  the  thermal  S2F10  conversion  mechanism 
discussed  in  the  previous  section,  it  was  necessary  to  investigate  the  pyrolysis  of 
S2F10  on  heated  metal  surfaces  at  low  pressure  and  at  different  temperatures.  The 
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Figure  20.  Single  ion  chromatogram  (ion  86)  of  4ppb  S2F10  in  200  kPa  SFe- 
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Figure  21.  Electron-impact  mass  spectra  of  1500 ppm  SFe  in  argon  buffer  gas  at  (a)  30°C 
and  (b)  220°C. 

specific  purposes  of  this  investigation  were  to:  1)  show  evidence  that  it  is  surface 
chemistry  that  is  involved  in  the  S2F10  conversion:  2)  indicate  the  temperature  range 
over  which  this  conversion  is  effective;  3)  demonstrate  that  SFe  does  not  chemically 
decompose  in  this  temperature  range:  and  4)  provide  an  independent  verification  that 
SOF2  is  indeed  one  of  the  products  of  the  S2F10  conversion  process. 

An  Extrel  quadrupole  mciss  spectrometer  was  used  for  these  measurements,  and  the 
sample  gas  was  admitted  to  the  MS  via  a heated  stainless  steel  bellows  connected  to  a 
variable  leak  valve.  The  temperature  of  the  bellows  could  be  varied  from  20  to  300°  C 
and  the  mass  spectra  could  be  observed  as  a function  of  the  temperature  of  the  inlet 
system.  The  gas  pressure  in  the  bellows  wa5  estimated  to  be  lower  than  1 millitorr 
(~0.13  Pa)  and  the  pressure  in  the  MS  ion-source  was  below  5 x 10“^  torr. 

Shown  respectively  in  figures  21  and  22  are  mass  spectra  obtained  from  mixtures  of 
1500 ppm  SFe  and  1500 ppm  S2F10  in  argon  for  different  temperatures  of  the  inlet 
bellows.  At  30  °C,  the  mass  spectra  of  S2F10  and  SFe  are  seen  to  be  nearly  identical 
in  agreement  with  previously  published  data  [30. 31. 37. 3S].  Other  than  minor 
differences  in  the  relative  intensities  of  the  SFJ  ions  (x  = l-4),  the  only  significant 
difference  between  the  mass  spectra  of  S2F10  and  SFq  is  the  appearance  of  a peak  at 
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Figure  22.  Electron- impact  mass  spectra  of  1500 ppm  S2F10  in  argon  buffer  gas  at  various 
temperatures. 
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ml z = 0^1  (SFj*^)  in  the  SFe  spectnim  which  is  not  present  in  S2F10  spectrum.  Peaks 
at  masses  48,  64,  67,  and  86  in  the  S2F10  spectrum  at  30  °C  are  due  to  small  impurities 
of  SO2  and  SOF2  present  in  the  S2F10  sample.  At  220^^0,  the  SFe  mass  spectrum 
remains  essentially  unchanged,  indicating  that  SFe  does  not  undergo  a measurable 
chemical  transformation  in  the  heated  inlet  tube.  The  S2F10  spectrum,  on  the  other 
hand,  exhibits  a significant  change  for  temperatures  above  about  100  °C.  In  going  from 
30  °C  to  200  °C,  the  ratio  of  the  SFj  (mjz  = 89)  to  SFj  [m/z  = 127)  peaks  in  the 
S2F10  mass  spectrum  increases  indicating  conversion  of  S2F10  into  SF4  or  SF5.  Peaks 
associated  with  SOF2  at  mjz  = 48,  67,  and  86  become  more  prominent  above  100 °C 
and  dominate  the  spectrum  at  270^^0,  indicating  that  S2F10  is  nearly  completely 
destroyed  at  this  temperature  with  SOF2  as  the  predominant  stable  by-product. 

Tests  performed  using  different  inlet  tubes  appear  to  indicate  that,  at  any  given 

temperature,  the  conversion  of  S2F10  into  SOF2  becomes  more  efficient  as  the  effective 

surface  axea  encountered  by  the  incoming  gas  is  increased,  e.g.,  a bellows  gives  a more 

complete  conversion  than  a straight  tube  of  the  same  length.  This  trend  is  consistent 

with  chemical  conversion  via  a surface  catalyzed  reaction.  There  are  indications  that 

some  surface  materials  may  be  more  effective  than  others  in  converting  S2F10  to 

SOF2,  however,  an  extensive  investigation  appeared  to  be  unwarranted  in  view  of  the 

• 

relatively  high  efficiency  achieved  with  commonly  used  stainless  steel.  Temperature 

and  presence  of  adsorbed  water  appear  to  be  more  important  than  either  surface 

configuration  or  composition. 

The  detailed  mechanisms  of  S2F10  conversion  to  SOF2  on  heated  surfaces  are  not 
presently  known.  It  is  speculated  that  the  process  involves  thermally  induced 
dissociation  of  S2F10  along  the  weak  sulfur-sulfur  bond  followed  by  a fast  reaction 
of  SF5  with  H2O  adsorbed  on  the  surface,  i.e., 

S2F10  — ^ SF5  + SF5 

SF5  + H2O  ->  SOF2  + 2HF  + F. 

Previous  work  [39-41]  has  shown  that  S2F10  is  thermally  unstable  in  the  gas  phase  at 
temperatures  above  100  °C.  The  thermal  unimolecular  gas-phase  decomposition  rate 
at  high  pressures  [39]  has  a temperature  ( T ) dependence  given  by: 

= 3 X 10^^exp(-46400/i?r),  (5) 

which  is  based  on  reaction  data  [39]  obtained  in  the  temperature  range  435-454.5  K 
(RT  is  given  in  units  of  calories  and  k in  s“^).  Assuming  this  expression  applies 
at  room  temperature,  the  half-life  of  S2F10  in  the  gas  phase  is  predicted  to  be 
about  10^  years.  The  half-life  drops  dramatically  to  10^  seconds  at  475  K.  However, 
the  S2F10  decomposition  rate  implied  by  the  present  mass-spectrometric  results  is 
significantly  higher  than  predicted  by  the  above  gas-phase  rate  and  unlikely  to  apply 
at  the  low  pressures  typically  encountered  here  in  the  heated  inlet  tube,  further 
supporting  the  likelihood  of  a surface  catalyzed  reaction  mechanism. 
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Figure  23.  Signal  intensity  of  several  selected  ions  from  a 1500  ppm  S2F10  sample  in  argon 
as  a function  of  time  as  the  inlet  temperature  is  rapidly  increased  from  room  temperature 
to  270°C. 


An  indirect  indication  of  the  role  of  adsorbed  water  was  obtained  from  measurements 
made  with  the  quadrupole  mass  spectrometer  under  conditions  where  the  inlet  line 
temperature  was  allowed  to  rise  relatively  rapidly  after  the  Ar/S2Fio  mixture  was 
introduced.  The  mtiss  spectra  obtained  in  this  way  differ  from  those  shown  in  figure  22 
which  correspond  to  equilibrated  conditions  achieved  after  keeping  the  inlet  line  at 
the  indicated  temperatures  for  many  minutes.  The  results  obtained  from  monitoring 
selected  ions  as  a function  of  time  under  the  '‘transient”  conditions  of  continuous 
rising  inlet-line  temperature  are  shown  in  figure  23.  Presented  in  this  figure  are  time- 
dependent  peak  intensities  for  ions  at  7^/2:  = S6  and  67  which  are  indicators  of  SOF2; 
mjz  ~ 127,  indicator  of  S2F10  and  SFe;  and  miz  = 89,  indicator  of  S2F10,  SFe,  and 
SF4.  Unlike  the  results  given  in  figure  22,  the  yield  of  SOF2  is  seen  to  go  through 
a maximum  at  a time  of  6 minutes,  which  corresponds  to  a temperature  of  160  °C. 
The  data  for  m/r  = 127  indicate  that  S2F10  is  nearly  completely  decomposed  at  this 
temperature,  i.e.,  the  intensity  for  this  ion  ceases  to  decrease  significantly  beyond 
this  temperature.  The  7n/z  = 127  signal  that  appears  at  higher  temperatures  is 
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undoubtedly  associated  with  SFe  which  can  be  a product  of  S2F10  decomposition.  As 
noted  by  Herron  [41]  (see  also  Ref.  [42]),  SF4  and  SFe  can  be  formed  in  the  gas  phase 
by  the  disproportionation  reaction 

SF5  + SF5  ^ SF4  + SFe. 

As  the  temperature  increases,  the  m/z  = 89  signal  does  not  fall  off  quite  as  rapidly 
as  the  mj z =■  127  signal,  thus  suggesting  formation  of  SF4  whose  primary  electron 
impact  ion  is  SFj . The  increase  in  SOF2  production  with  time  can  be  attributed  to  an 
increase  in  the  decomposition  rate  of  S2F10  and  an  increase  in  the  surface  hydrolysis 
rate  of  the  decomposition  products  as  the  surface  temperature  is  increased.  However, 
as  the  amount  of  water  adsorbed  on  the  surface  decreases  due  to  heating,  the  rate  of 
S2F10  conversion  into  SOF2  drops.  It  should  be  noted  that  at  the  highest  temperatures 
the  production  of  SOF2  has  decreased  with  no  corresponding  increase  in  any  of  the 
other  compounds  which  are  being  monitored.  This  suggests  that  at  high  temperatures 
the  S2F10  or  its  by-products  are  being  converted  into  compounds  which  are  presently 
unidentified. 

Although  there  is  evidence  [43]  that  SOF2  hydrolyzes  slowly  in  the  gas  phase  at  300  K 
and  may  react  with  H2O  on  surfaces,  little  is  known  about  its  decomposition  at  higher 
temperatures.  These  results  indicate  that  there  is  an  optimum  temperature  range  at 
which  conversion  of  S2F10  into  SOF2  is  most  efficient.  If  the  temperature  is  too  low, 
the  conversion  rate  is  low,  and  if  the  temperature  is  too  high,  SOF2  itself  may  be 
destroyed. 


2.2.5  Gas-phase  Impurities 

It  should  be  realized  that  if  gaseous  species  are  present  in  the  unknown  sample  that 
could  react  with  S2F10  or  its  decomposition  products  at  temperatures  above  150° C, 
uncertainties  can  be  introduced  in  the  quantitative  analysis  for  S2F10  by  the  method 
proposed  here  if  the  reference  sample  does  not  also  contain  these  species  in  nearly 
the  same  amounts.  It  is  known  [44-47]  that  S2F10  can  react  at  high  temperatures 
with  numerous  gases  such  as  NO,  NO2,  NO2CI,  (CF3)2NONO,  CH3ONO,  Fe(CO)5, 
ICl,  CeHg,  C2H4,  (CH3)2NH,  BCI3,  and  NH3.  If  it  is  suspected  that  gases  may  be 
present  in  the  unknown  sample  that  can  react  with  S2F10  at  temperatures  above 
150° C,  then  care  must  be  taken  to  assess  the  influence  of  these  gases  on  the  GC-MS 
response  to  S2F10.  Ideally  the  reference  sample  should  have  the  same,  or  nearly  the 
same  composition  as  the  unknown  sample.  It  was  found  here  that  the  response  of  the 
GC-MS  to  S2F  10  was  not  significantly  influenced  by  the  presence  of  the  oxyfluorides 
SOF2,  SO2F2,  and  SOF4  which  are  commonly  formed  together  with  S2F10  during 
decomposition  or  oxidation  of  SFe. 

Of  particular  importance  is  the  presence  of  water  vapor  as  a contaminant  in  a sample. 
The  gas-phase  water  vapor  present  in  the  sample  can  affect  the  amount  of  adsorbed 
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water  on  the  inlet  tube  surface  at  any  given  time  and  thus  the  S2F10  to  SOF2 
conversion  efficiency.  Low  water-vapor  concentrations  in  gas  samples  did  not  show 
significant  effects.  However,  it  was  observed  that  the  addition  of  excessive  amounts 
of  water  vapor  to  an  S2F10  sample  can  noticeably  decrease  the  amount  of  SOF2 
detected  from  S2F10  conversion  at  surface  temperatures  above  about  120° C.  This  is 
not  understood  but  may  be  due  to  mechanisms  which  either  inhibit  the  conversion  if 
gas-pha^e  water  is  present,  or  increase  the  rate  of  SOF2  decomposition  (for  example 
by  hydrolysis  [35])  at  elevated  temperatures.  The  reactivities  of  gas-phase  S2F10  and 
SOF2  toward  water  vapor  at  temperatures  above  100°  C are  presently  unknown  and 
require  more  investigation. 


2.2.6  Reference  Sample  Stability 

Surface  conversion  of  S2F10  into  SOF2  also  occurs  at  room  temperature,  although 
at  a much  slower  rate  than  observed  at  elevated  temperatures.  This  decomposition 
raises  questions  about  the  long-term  reliability  of  reference  gas  samples  containing 
known  amounts  of  S2F10  that  are  needed  for  calibration  of  the  GC-MS.  In  order  to 
assess  the  stability  of  S2F10  in  reference  samples,  a series  of  tests  were  performed  to 
monitor  the  content  of  S2Fio-containing  cylinders  over  relatively  long  periods  of  time. 
Preliminary  reports  of  these  tests  have  already  been  given  [10,24]. 

It  was  found  that  the  rate  of  S2F10  decomposition  depended  on  such  factors 
as  temperature,  cylinder  material  and  surface-to-volume  ratio,  and  moisture 
content.  Data  indicating  the  temperature  and  moisture  content  dependencies  of 
S2F10  decomposition  in  150  ml  stainless-steel  cylinders  are  shown  respectively  in 
figures  24  and 25.  Shown  in  these  figures  are  measured  S2F10  concentrations  in  argon 
as  a function  of  time  after  mixture  preparation.  The  measurements  were  performed 
using  the  GC-MS  method  described  above,  and  the  “aged”  samples  were  compared 
with  reference  samples  prepared  immediately  prior  to  the  measurements.  The  “wet” 
sample  (H2O  added)  results  in  figure  25  were  obtained  from  a sample  prepared  by 
injecting  1 pi  of  liquid  water  into  the  cylinder  under  vacuum  prior  to  introducing  the 
Ar/S2Fio  gas  mixture.  The  data  presented  in  figures  24  and  25  clearly  demonstrate 
that  the  S2F10  decomposition  rate  increases  with  increasing  temperature  and  moisture 
content. 

The  products  of  the  S2F10  decomposition  were  found  to  include  SOF2  and  SFe,  and 
thus  it  is  speculated  that  the  decomposition  process  involves  a surface  catalyzed 
reaction  with  adsorbed  H2O  that  may  be  similar  to  that  observed  at  higher 
temperatures  in  the  heated  MS  gas  inlet  tube.  The  SOF2  produced  from  S2F10  also 
appears  to  react  slowly  on  the  stainless-steel  surface  and  convert  to  other  products 
that  have  yet  to  be  identified.  Once  S2F10  disappears  from  the  gas,  SOF2  will  also 
eventually  disappear  [38].  From  analyses  of  the  decomposed  gas  using  MS  and  GC- 
MS,  it  has  not  been  possible  to  account  for  all  of  the  sulfur  originally  present  in  the 
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ELAPSED  TIME  (hours) 


Figure  24.  Low  temperature  dependence  for  decomposition  of  500  ppm  S2F10  samples  in 
argon  stored  in  150  ml  stainless  steel  cylinders. 


ELAPSED  TIME  (days) 


Figure  25.  Decomposition  of  S2F10  as  a function  of  time  for  two  1400  ppm,  S2F10  samples 
in  1 atm  Ar  in  150  ml  stainless  steel  cylinders.  For  one  sample  (□),  1 /il  of  liquid  water  was 
injected  into  the  cylinder  prior  to  filling  with  gas. 


2.3  Stochastic  Properties  of  Pulsating  Partial  Discharges 


33 


SoFio-  The  only  identified  gaseous  compound  resulting  from  S2F10  decomposition 
that  contains  sulfur  and  exhibits  long-term  stability  is  SFe-  Generally,  less  than  15  % 
of  the  sulfur  can  be  accounted  for  by  this  product. 

On  the  basis  of  the  tests  that  were  conducted,  long-term  storage  of  S2F10  reference 
samples  in  the  gas  phase  cannot  be  recommended.  It  was  found  that  samples 
stored  under  even  the  most  favorable  conditions  exhibited  significant  decomposition 
over  sufficiently  long  periods  of  time.  It  is  recommended  that  reference  samples 
be  prepared  immediately  prior  to  gas  analysis,  under  conditions  that  are  as  dry  as 
possible,  and  in  cylinders  that  are  as  large  as  possible. 

Preliminary  tests  of  commercial  grade  SFe  (99.8%  purity)  have  indicated  that  S2F10 
does  not  readily  decompose  in  SFs-filled  cylinders  when  compressed  into  the  liquid 
phase.  Analysis  showed  that  S2F10  was  present  as  an  impurity  at  a level  of  about 
40ppb  in  a 12-year  old  sample.  It  should  be  noted  that  S2F10  is  a possible  by-product 
from  the  commercial  production  of  SFe  and  can  be  present  as  a trace  liquid  impurity 
in  the  liquified  SFe  which  exists  in  pressurized  cylinders.  The  results  of  analyzing 
SFe  stored  as  a liquid  suggest  that  significant  improvements  might  be  achieved  in 
maintaining  sample  stability  by  compressing  reference  samples  into  the  liquid  phase. 


2.3  Stochastic  Properties  of  Pulsating  Partial  Discharges 

2.3.1  Basic  Considerations 

Partial-discharge  (P-D)  measurements  have  long  been  recognized  as  an  important 
part  of  quality  control  for  high-voltage  apparatus.  There  is  also  increasing  interest  in 
the  use  of  P-D  measurement  as  a diagnostic  tool  for  assessing  insulation  performance 
and  integrity  when  subjected  to  electrical  and  mechanical  stress.  The  measurement 
of  P-D’s  has  been  applied  to  a wide  range  of  electrical  apparatus  including  rotating 
machines,  power  transformers,  and  both  low-  and  high-voltage  cables  [48]. 

The  P-D  phenomenon  under  consideration  here  is  a self-quenching  discharge  that  is 
localized  in  regions  of  high  electrical  stress  within  insulating  systems.  The  high-stress 
regions  might,  for  example,  be  associated  with  voids  or  “trees*’  in  solid  dielectrics, 
metal  particles  or  sharp  points  on  either  solid  conductors  or  insulators,  or  the  presence 
of  interfaces  between  different  types  of  insulating  or  conducting  media.  Partial 
discharges  are  self-quenching  because  the  surface  or  space  charge  produced  during 
the  discharge  activity  generally  accumulates  and  acts  to  reduce  the  local  electric-field 
magnitude  to  a level  below  which  the  discharge  can  be  sustained. 

Because  of  its  self- quenching  nature,  the  P-D  phenomena  is  inherently  pulsating  and 
will  manifest  itself  as  current  pulses  in  the  external  circuit.  Partial-discharge  activity 
occurs  in  regions  that  are  at  least  partially  in  the  gas  phase  such  as  occlusions  in  solids 
or  “bubbles”  formed  by  vaporization  of  a liquid.  The  partial-discharge  phenomena 
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that  occur  around  conductors  in  air  or  other  gases  is  often  referred  to  as  corona. 
It  should  be  noted  that  corona  discharges  can,  under  some  conditions,  exhibit  a 
non-pulsating  characteristic  for  dc  voltages  which  is  associated  with  the  formation 
of  a steady  glow  region  [49,  50].  Because  this  type  of  corona  discharge  is  steady  and 
therefore  not  self- quenching,  it  does  not  fit  into  the  category  of  P-D  phenomena  under 
consideration  here. 

Partial  discharges  can  be  observed  by  electrical,  optical,  or  acoustical  measurements 
that  detect  and  record  pulses.  The  determination  of  P-D  pulse  characteristics 
has  been  the  subject  of  numerous  previous  investigations  [51].  Techniques  have 
been  developed  to  measure  such  parameters  as  mean  P-D  pulse  amplitude  and 
repetition  rate  [52],  pulse-height  distribution  [53-55],  pulse  shapes  [56],  correlation  of 
pulse  amplitude  with  phase  of  an  ac  voltage  [57,  59],  and  correlation  between  optical, 
electrical,  or  acoustical  pulses.  Because  of  the  complex,  stochastic  nature  of  P- 
D phenomena,  the  results  of  such  measurements  have  often  been  difficult,  if  not 
impossible,  to  interpret  in  terms  of  meaningful  physical  models.  The  measurement 
scheme  discussed  here  is  an  extension  of  a technique  previously  used  to  provide  more 
refined  information  about  the  statistical  characteristics  of  P-D  pulses.  From  the 
additional  information  acquired  by  this  method,  it  is  possible  to  gain  more  insight 
into  the  physical  bases  of  the  phenomena  which  might  prove  useful  in  the  design  of 
P-D  measurement  systems  of  the  future. 

Partial  discharge  phenomena  fit  into  a category  of  stochastic  processes  for  which 
memory  effects  play  an  important  role.  The  existence  of  memory  is  an  inherent 
property  of  P-D  phenomena  because  the  occurrence  of  a discharge  pulse  is  associated 
with  the  formation  of  residual  ion  space-charge,  surface  charges,  and  excited 
metastable  species  which  can  affect  the  initiation  and  growth  of  subsequent  discharge 
pulses.  As  mentioned  above,  the  pulsating  characteristic  of  P-D’s  in  various  types  of 
insulating  media  can  be  attributed  to  effects  of  space  or  surface  charges.  When  the 
accumulation  of  discharge-generated  space  or  surface  charge  causes  the  field  to  fall 
below  a level  sufficient  to  sustain  the  discharge,  it  ceases.  Another  P-D  pulse  will 
not  occur  at  this  site  until  the  local  field  is  restored  such  as  can  result  from  charge 
migration  or  changes  in  applied  voltage. 

Because  of  memory  effects,  correlations  may  exist  among  successive  discharge  pulses 
and  measurement  of  these  correlations  can  provide  useful  information  on  the  nature 
of  P-D  phenomena.  Recognition  of  this  fact  motivated  the  development  of  a system 
to  extract  information  about  memory  effects.  This  system  is  described  in  a recent 
publication  [60]  and  has  been  applied  to  an  investigation  of  the  stochastic  behavior 
of  Trichel-pulse,  negative-corona  discharges  in  electronegative  gases  [52,  61, 62]  and 
dielectric-barrier-type  discharges  in  air  which  will  be  discussed  in  the  next  three 
sections  [13]. 

It  is  convenient  to  represent  the  P-D  phenomenon  as  a random  point  process  [63]  in 
which  each  pulse  is  identified  by  a time  and  is  “marked”  with  a particular  property, 
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REPRESENTATION  OF  A MARKED 
RANDOM  POINT  PROCESS 
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Figure  26.  Diagrammatic  representation  of  the  marked  random  point  process 
corresponding  to  the  Trichel-pulse  discharge  phenomenon.  The  are  pulse  time 

separations  and  the  are  pulse  amplitudes. 

in  this  case  the  pulse  amplitude  as  indicated  in  figure  26  The  stochastic  process  is 
thus  defined  by  the  set  of  ordered  pairs  ^n},  n = 1,  2,  3, . . . Because  of  memory 
effects,  qn  and  tn  can  depend  on  the  set  for  all  j < n.  The  process  can  also 

be  defined  by  the  set  {^i,  n = 2,3,...,  j,  where  Af„_i  = — tn-i, 

which  is  more  appropriate  to  this  discussion  since  it  is  the  time  separations  between 
pulses  that  are  either  measured  or  specified  by  the  nieasurement  process.  The  pulse 
amplitude  was  chosen  as  the  relevant  mark  because  it  is  a measure  of  discharge 
intensity.  In  defining  the  point  process,  one  could  also  consider  other  marks  such  as 
the  pulse-shape  parameters. 

If  memory  is  important  in  the  P-D  process,  the  variables  ^n,Afn,g„_i  and  Afn-i 
associated  with  adjacent  events  may  not  be  independent.  For  example,  the  amplitude, 
may  depend  on  the  time  separation,  Af„_i.  Indeed,  this  is  found  to  be  the  case  for 
Trichel-pulse  discharges  [12,13,60-62,64].  The  dependence  of  Af„  on  Qn  or  Afn_i  is 
an  indication  of  the  influence  of  residuals  from  the  previous  pulse  in  either  enhancing 
or  retarding  the  initiation  of  the  subsequent  pulse.  Likewise,  the  dependence  on 
Afn-i  or  is  an  indication  of  the  influence  of  residuals  from  the  previous  pulse  on 
the  growth  of  the  subsequent  pulse. 

Because  of  these  dependencies,  the  unconditional  pulse-amplitude  and  pulse  time- 
separation  distributions  (po(qn)  arid  po(Afn))  are  not  independent.  If,  for  example, 
an  external,  radiation  source  is  used  to  help  initiate  P-D’s,  po(Atn)  will  change  with 
the  intensity  of  this  source  since  the  probability  of  discharge  initiation  changes  with 
intensity.  The  changes  in  po(Af„)  will  then  be  reflected  in  po(qn)  if  qn  depends  on 
Af„_i.  Results  from  measurements  of  only  one  distribution,  po(qn)  or  po(Atn),  as 
are  obtained  with  most  conventional  P-D-measurement  systems,  will  be  difficult  to 
interpret  if  memory  effects  are  significant. 

The  ability  to  detect  correlations  among  successive  discharge  pulses  is  a valuable 
diagnostic  that  could  be  incorporated  into  existing  P-D  measurement  systems 
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BLOCK  DIAGRAM  OF  THE  STOCHASTIC  MEASUREMENT  SYSTEM 


Figure  27.  Diagrammatic  representation  of  system  for  measuring  pulse-amplitude  and 
time-separation  distributions.  In  this  diagram,  the  relevant  components  are  identified  as 
follows:  MCA  - Multichannel,  analyzer;  TAG  - time-to-amphtude  converter;  SCA  - single 
channel  analyzer;  DDG  - digital-delay  generator;  A - amplifier;  G - gate;  and  S - switch. 

either  directly  through  filtering  techniques  described  here,  or  indirectly  through 
off-line  data  analysis.  The  NIST  system  allows  “real-time”  measurement  of  a 
set  of  conditional  pulse-amplitude  and  pulse-time-separation  distributions,  namely 
Pi(^nkn-i);  j > 1;  pi(At„|(f),  ^ or  At„_i;  and  p2(gn|Atn-i,  0^ 

(f  = qn-i  or  Atn-2-  These  distributions  are  defined,  for  example,  such  that 
pi{qn\Atn-i)dqn  is  the  probability  that  the  nth  P-D  pulse  has  an  amplitude  in  the 
range  qn  to  + dqn  if  it  is  separated  from  the  previous  pulse  by  a fixed  time  Atn_i, 
and  p2(5n|Af„_i,  is  the  same  but  with  both  Afn-i  and  qn-i  fixed. 

By  comparison  of  the  measured  conditional  and  unconditional  pulse-height  and 
time-separation  distributions  it  is  possible  to  infer  immediately  whether  or  not  two 
variables  are  correlated  and  therefore  dependent.  If  it  is  found,  for  example,  that 
Pi{qn\Atn-i)  ^ Poiqn)  at  least  some  allowed  values  of  qn  and  for  all  allowed  Atn_i 
such  that  po{Atn-i)  i=-  0,  then  qn  and  Atn_i  are  dependent  variables.  Similarly,  if 
P2(^n|Af„_i, ^„_i)  ^ pi(qn\Atn-i)  for  cillowed  g„_i,  then  a dependence  must  exist 
between  and  qn-i  at  a fixed  At„_i.  One  can  also  use  the  measured  conditional 
distributions  to  compute  correlation  coefficients  that  provide  a quantitative  measure 
of  the  extent  to  which  any  two  variables  in  the  set  {gi,  g„,  Atri_i  }j  are  correlated. 

Measurement  of  meaningful  conditional  distributions  require  the  P-D  process  to 
be  stationary,  or  at  least  quasistationary.  If  external  factors  that  affect  P-D 
behavior,  such  as  the  size  of  a void,  or  the  electron  emission  properties  of  a surface, 
change  during  the  measurement,  the  data  obtained  for  conditional  distributions 
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could  give  false  indications  of  possible  correlations.  Tbe  extent  to  which  external 
factors  are  important  is  perhaps  best  revealed  by  measurements  of  the  unconditional 
distributions  po(?n)  and  po{Atn)  which  can  be  performed  relatively  more  quickly  than 
the  measurements  of  conditional  distributions.  The  Trichel-pulse  corona  discharge 
considered  below  is  an  example  of  a P-D  phenomenon  that  is  sufficiently  stationary 
so  as  to  allow  measurement  of  all  the  conditional  distributions  for  which  the  system 
is  designed. 

2.3.2  Measurement  System 

The  modifications  of  the  Stochastic  Analyzer  for  Pulsating  Phenomena  (SAPP) 
originally  described  by  Van  Brunt  and  Kulkarni  [60]  have  been  made  to:  1)  allow 
direct  real-time  measurements  of  the  conditional  pulse- amplitude  distributions 
Pi{qn\Atn-j)^  j > 2 andp2(^niA  t„_i,  A t„_2);  2)  eliminate  intermediate  pulse  errors 
associated  with  measurement  of  the  p2  distributions;  and  3)  reduce  background 
associated  with  accidental  sampling  of  the  tails  of  discharge  pulses.  A diagram  of 
the  measurement  system  is  shown  in  figure  27.  It  should  be  noted  that  another 
version  of  this  system  has  recently  been  constructed  at  NIST  which  can  be  applied 
to  investigation  of  P-D’s  generated  by  ac  voltages.  This  new  system,  which  will  be 
described  later,  allows  measurement  of  distributions  that  are  conditioned  on  a set  of 
selectable  phase  intervals  of  the  applied  voltage. 

The  system  in  figure  27  may  be  configured  to  measure  any  of  the  indicated 
distributions  by  appropriate  positioning  of  the  switches  (S1-S5).  Data  from  all 
measurement  configurations  are  accumulated  by  a 256-channel,  computer-controlled 
multi-channel  analyzer  (MCA). 

The  system  was  first  applied  to  an  investigation  of  the  stochastic  behavior  of  negative 
corona  (Trichel  pulses)  generated  using  a point-plane  electrode  gap  configuration  in 
electronegative  gas  mixtures  [49].  The  test  gap  is  indicated  in  figure  27.  The  corona- 
discharge  current  pulses  are  detected  electrically  using  a preamplifier  (Al)  connected 
to  an  impedance  Z that  is  in  series  with  the  discharge  gap.  The  detected  pulses 
are  then  sent  to  a variable  gain  amplifier  (A2)  before  being  routed  to  the  appropriate 
circuit  path.  For  measurement  of  po(^n),  the  discharge  pulses  from  A2  are  fed  directly 
to  the  MCA  through  gate  G3  which  is  kept  open  by  proper  positioning  of  S4. 

The  distribution  po(Atn)  is  measured  by  sending  the  pulses  from  A2  to  a pulse  sorter 
used  to  trigger  two  time-to-amplitude  converters  (TACl  and  TAC2)  which  generate 
output  pulses  of  amplitude  proportional  to  the  time  between  P-D  pulses.  The  outputs 
of  TACl  and  TAC2  are  fed  to  the  MCA  through  gate  G2  by  appropriate  positioning  of 
S3.  As  shown  in  ref.  [60],  this  arrangement  permits  the  recording  of  every  sequential 
time  interval,  provided  all  time  separations  exceed  the  TAG  reset  time  of  50  ps. 
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Tlie  measurement  of  pi(g„|Afn_i)  requires  use  of  a logic-controlled  digital-delay 
generator  (DDG)  to  restrict  the  transfer  of  pulses  from  A2  to  the  MCA  through  G3  to 
a narrow  preselected  time  interval  Atn-i  ±^(Af„_i)/2  for  which  Atn-i  » (5(Af„_i). 
To  measure  p2(^n|Af„_i,  gn_i),  a single-channel  analyzer  (SCAl)  is  inserted  between 
A2  and  the  input  to  the  DDG  time-interval  control  logic  circuit  by  proper  positionings 
of  Si  and  S5.  The  SCAl  restricts  triggering  of  the  DDG  to  pulses  having  amplitudes 
within  the  selected  range  ±S{qn-i)/2  for  which  g„_i  >>  S(qn-i)  [60].  To  insure, 
in  the  case  of  p2  measurements,  that  the  time-interval  control  logic  properly  inhibits 
the  opening  of  G3  if  intermediate  pulses  occur  before  the  selected  time  window,  it 
is  necessary  to  modify  the  circuit  previously  described  [60]  by  adding  the  amplifier 
A3  with  input  at  /.  This  amplifier  is  identical  to  that  defined  by  transistor  Ti  in 
figure  16  of  ref.  [60].  In  the  present  configuration,  the  input  / is  connected  directly 
to  the  delay  Ti  of  figure  16  (ref.  [60])  which  has,  in  turn,  been  disconnected  from  Ti. 

Unlike  the  system  previously  described  [60],  the  present  system  allows  direct 
measurement  of  p2(^n|A^n-i?  For  measurement  of  this  distribution,  the 

s^vitches  Si  and  S5  are  positioned  so  that  the  input  to  SCAl  is  derived  from  TACl. 
The  time  interval  measured  by  TACl  is  restricted  by  SCAl  to  lie  within  a narrow 
range  Af„_2  ± (5(At„_2)/2  for  triggering  the  DDG,  which,  in  turn,  defines  At„-i* 

By  the  indicated  positionings  of  switches  Si,  S2,  and  S5,  it  is  possible  to  measure  the 
conditional  pulse  time-separation  distributions  pi(Af„|Atn_i)  and  pi{Atn\qn)-  For 
measurement  of  pi(Af„lAf„_i),  the  output  of  TACl  is  fed  to  the -input  of  the  single- 
channel analyzer  SCA2  which  thus  defines  Af„_i  and  triggers  TAC3  for  measurement 
of  the  next  time  interval  Af„.  In  the  case  of  the  pi(Af„|^ri)  measurement,  SCA2  is 
connected  directly  to  amplifier  A2  thereby  used  to  define  qn. 

The  present  system  also  differs  from  that  previously  described  in  that  the  pi(^n|Ai„-j) 
logic  circuit  now  contains  a pulse  counter  that  can  be  set  for  any  j > 2.  By  adjustment 
of  the  delays  inherent  to  the  time-interval  control  logic  and  with  minor  adjustments  in 
the  shape  of  the  gate  pulse  from  this  circuit,  it  has  been  possible  to  largely  eliminate 
the  background  problem  illustrated  by  figure  19  of  ref.  [60]. 

In  the  present  measurement  system,  the  pulse  detection-circuit  conditions  are  similar 
to  those  described  in  ref.  [60]  so  that  the  measured  discharge-pulse  amplitude  is 
proportional  to  the  net  charge  transported  in  the  discharge  given  by 


(6) 


where  is  the  instantaneous  discharge  current  at  time  t for  the  7ith  pulse.  Given 
the  impulse  response  h(t  — t')  of  the  detection  circuit  (defined  here  mainly  by  the 
filter  characteristics  of  Z and  Al),  the  observed  pulse  signal  is  given  by 


it)  = /_ 


— oo 


in  {t')  h{t  — t')  dt' . 
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The  impulse  response  in  the  present  measurement  system  has  a width,  u;,-,  of  about 
1.5  ps  whereas  the  intrinsic  width,  Wp^  of  a typical  Trichel  pulse  is  known  [56]  to  be 
on  the  order  10  to  50  ns,  Therefore  the  condition  Wi  » Wp  allows  the  approximation 

i„(<)  ~ - t„),  (S) 

where  8(t  — f„)  is  the  Dirac  delta  function.  Using  this  in  equation  (7)  gives 

9n(<)  = (9) 

This  means  that  the  shapes  of  the  observed  Trichel  pulses  in  the  present  experiment 
are  governed  primarily  by  the  shape  of  h{t  — and  thus  the  maximum  of 
denoted  here  by  is  directly  proportional  to  Q'^.  Therefore,  can  be  expressed  in 
units  of  charge,  and  in  the  present  case,  it  is  convenient  to  use  picocoulombs  (pC). 
The  method  for  calibrating  the  pulse-height  measurements  in  terms  of  charge  is  the 
same  as  described  previously  [55]. 

2.3.3  Examples  of  Results  from  Investigation  of  Trichel  Pulses 

The  Trichel-pulse  discharge  [65]  was  selected  for  initial  test  of  the  SAPP  measurement 
concept  because  its  properties  are  reasonably  well  understood  [66-70]  and  it  can 
be  easily  generated  and  controlled  in  a simple  point-plane  gap.  With  some  care, 
the  Trichel-pulse  discharge  can  be  maintained  in  a stable  condition  for  relatively 
long  periods  of  time.  The  phenomenon  is  shown  [12]  to  be  a clear  example  of 
a non- Markovian,  marked  random  point  process  in  which  memory  effects  play  an 
important  role.  Strong  correlations  are  found  to  exist  among  the  amplitudes  and 
time  separations  of  successive  discharge  pulses  which  can  be  interpreted  in  terms  of 
the  influences  of  ion  space  charge  and  metastable  species  from  previous  pulses  on  the 
initiation  and  growth  of  subsequent  pulses. 

As  an  example  of  the  richness  of  information  that  can  be  extracted  with  the  present 
measurement  system,  results  obtained  for  self-sustained  Trichel  pulses  in  a neon- 
oxygen  gas  mixture  are  shown  in  figures  28-32.  Data  for  the  different  sets  of 
distributions  were  obtained  at  different  times  so  that  the  cathode  surface  conditions 
that  apply,  for  example,  to  the  results  in  figure  28  differ  slightly  from  those  that  apply 
to  figure  29.  This  accounts  for  the  difference  in  po(qn)  shown  in  these  two  figures. 
A detailed  interpretation  of  the  results  presented  here  would  go  beyond  the  scope  of 
this  report.  However,  salient  features  of  the  data  and  certain  important  conclusions 
that  can  be  derived  therefrom  should  be  noted.  A more  complete  discussion  of  the 
physical  bases  for  the  stochastic  behavior  of  the  Trichel-pulse  phenomenon  has  been 
given  by  Van  Brunt  and  Kulkarni  [12]. 

The  fact  that  the  second-order  conditional  distributions,  p2,  differ  from  the 
corresponding  first-order  distributions,  pi,  which  in  turn  differ  from  the  corresponding 
unconditional  distributions,  po,  indicates  unequivocally  that  the  set  of  variables 
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CONDITIONAL  AND  UNCONDITIONAL  PULSE-AMPLITUDE 

DISTRIBUTIONS 


qn(pC) 


Figure  28.  Measured  unconditional  and  conditional  discharge  pulse- amplitude 
distributions  po(?n),  Pi(?n| )?  and  p2(?n|At„_i, g„_i)  for  the  indicated  fixed  values 
for  At„_i  and  Qn-i-  The  distributions  have  been  normalized  to  the  maxima.  These  results 
are  for  negative-corona  pulses  in  a neon-oxygen/gas  mixture. 
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Figure  29.  Measured  unconditional  and  conditional  discharge  pulse- amplitude 
distributions  po(?n),  Pi{Qn\^tn-i)-,  and  P2(5nl At„_2)  for  the  indicated  fixed  values 
for  Atn-i  and  Atn-2-  The  distributions  have  been  normalized  to  the  maxima.  These  results 
are  for  negative-corona  pulses  in  a neon-oxygen/gas  mixture. 
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CONDITIONAL  AND  UNCONDITIONAL  PULSE 
TIME-SEPARATION  DISTRIBUTIONS 


Atn(|XS) 

Figure  30.  Measured  unconditional  and  conditional  discharge  pulse-time-separation 
distributions  po{Atn),  pi(At„|g„),  and  At„_i)  for  the  indicated  fixed  values  for  Qn 

and  At„_i.  The  distributions  have  been  normalized  to  the  areas  under  the  curves.  These 
results  are  for  negative-corona  pulses  in  a neon- oxygen /gas  mixture. 
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Figure  31.  Measured  unconditional  and  conditional  discharge  pulse-amphtude 
distributions  Pi(5n|^^n-2)5  Pi(?n|?n-i)  for  the  indicated  hxed  values  for  At„-2 

and  Qn-i-  The  distributions  have  been  normalized  to  the  areas  under  the  curves.  These 
results  are  for  negative-corona  pulses  in  a neon- oxygen /gas  mixture. 
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Figure  32.  Measured  unconditional  and  conditional  discharge  pulse- amplitude 
distributions  'Po(qn)’,  Piiqn\^in-3)-,  and  pi(g„|  At„_4)  for  the  indicated  fixed  values  for  Atns 
and  Atn-4-  The  distributions  have  been  normalized  to  the  maxima.  These  results  are  for 
negative-corona  pulses  in  a neon- oxygen /gas  mixture. 
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{A^n,  Atn_i,  Qn-i,  At„_2, ' * '}  associated  with  adjacent  pulses  are  not  independent. 
For  example,  it  is  seen  from  figures  28  and  29  that  qn  has  a strong  positive  dependence 
on  A^n-i*  This  dependence  can  be  related  to  the  influence  of  negative-ion  space 
charge  from  previous  pulses  in  suppressing  the  magnitude  of  the  electric  field  at 
the  cathode  when  the  next  pulse  develops.  It  is  also  seen  from  figure  28  that  the 
amplitude,  g„,  of  a pulse  can  be  either  positively  or  negatively  dependent  on  the 
amplitude  of  the  previous  pulse.  The  sign  of  this  dependence  can  be  explained  in 
terms  of  the  competing  effects  of  negative-ion  space  charge  and  metastable  species 
in  respectively  retarding  or  enhancing  the  growth  of  the  next  pulse.  The  negative 
dependence  of  Atn  on  qn  (and  Af„_i)  implied  by  the  conditional  time-separation 
distributions  shown  in  figure  30  can  be  understood  in  terms  of  the  influence  of 
metastable  species  from  the  previous  pulse  in  enhancing  the  probability  for  initiating 
the  next  pulse  by  ejecting  electrons  from  the  cathode  surface  during  field- assisted 
quenching. 

Because  of  the  correlations  among  the  amplitudes  and  time  separations  of  successive 
pulses,  the  distributions  shown  in  figures  28-30  are  all  related.  It  can  be  shown, 
for  example,  from  the  law  of  probabilities  that  poiqn),  Po(Atn),  and  pi Atn-i ) are 
related  by  the  integral  expression 

rod 

Poiqn)  = / po(A^n_i)pi(^„|Af„_i)d(A^„_i),  (10) 

Jo 

and  the  distributions  po(gn),Po(At„),pi(5„|Afn-i),Pi(Af„|g„),  andp2(^n|<?n,  At„_i) 
are  related  by 

roo 

Po(Atn-l)pi  (g„|Af„_i)  = / Po{qn-l)Pl  (Atn-l\qn-l)  P2  ( kn-1 , At  „_i  ) . 

Jo 

(11) 

Equation  (10)  indicates  that  if  qn  is  dependent  on  Afn-i,  then  any  externally- 
induced  change  in  the  time-interval  distribution,  po(Afn),  will  necessarily  change 
the  amplitude  distribution,  po(?n)- 

Since,  as  seen  from  the  data  in  figures  29, 31,  and,  32,  the  profiles  for 
P2(?n|Af„_i,  Afn_2)  and  for  pi(^„|Af„_j),  j = 2,3,4  do  not  match  the  profile  for 
Po(qn),  it  can  be  stated  that  depends  on  Afn-j,i  > T and  therefore,  the  process 
is  one  for  which  memory  extends  back  in  time  beyond  the  most  recent  event,  i.e., 
the  process  is  non-Markovian.  This  observation  is  consistent  with  results  reported 
in  the  recent  work  of  Steiner  [64].  It  can,  in  fact,  be  shown  [12]  that  because  of  the 
relatively  strong  dependence  of  on  both  A^n-i  and  ^n-i,  it  is  possible  for  memory 
to  propagate  indefinitely  back  in  time. 
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Figure  33.  Electrode  configuration. 

2.3.4  Influence  of  a Dielectric  Barrier  on  the  Stochastic  Behavior  of  Trichel-Pulse 
Corona 

The  investigation  discussed  here  was  undertaken  to  examine  the  influence  that  a thin 
solid  polytetrafluoroethylene  (PTFE)  dielectric  barrier  has  on  the  stochastic  behavior 
of  negative  (Trichel)  pulse  corona  discharges  in  air.  Although  the  effects  of  solid 
dielectrics  on  corona  pulses  have  been  examined  in  previous  experiments  [71],  there 
were  no  attempts  to  characterize  the  stochastic  behavior  of  the  phenomenon  in  these 
investigations.  It  is  shown  here  from  the  measurements  of  conditional  pulse-amplitude 
distributions  that  corona-induced  charging  of  dielectric  barriers  placed  on  the  planar 
electrode  introduces  another  memory  effect  that  becomes  increasingly  important  as 
the  point-to-plane  gap  spacing  is  reduced.  It  is  also  shown  that,  because  the  PTFE 
dielectric  acquires  a “quasi-permanent”  charge  from  the  corona,  there  is  a critical  gap 
spacing  below  which  the  Trichel- pulses  appear  to  cease. 

The  experiments  were  carried  out  using  the  SAPP  technique  described  above.  The 
electrode-gap  configuration  used  to  obtain  the  results  presented  here  is  shown  in 
figure  33.  A PTFE  dielectric  of  1.0  mm  thickness  and  of  diameter  D was  placed  on  the 
plane  electrode  and  centered  on  the  point-to-plane  axis.  The  point-to-plane  spacing, 
d,  could  be  varied  from  0 to  10  cm.  A very  sharp  stainless-steel  point  electrode  with 
a radius  of  curvature  at  the  tip  of  less  than  0.01  mm  was  used  to  insure  that,  at  a 
sufficiently  high  gap  voltage,  Idi,  the  discharge  could  be  sustained  by  field-induced 
electron  emission.  The  measurements  described  here  were  performed  using  “room 
air”  in  the  gap.  The  corona-discharge  current  pulses  were  detected  electrically  using 
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Figure  34.  The  measured  critical  point-to-plane  gap  spacing,  below  which  the 
Trichel-pulse  discharge  ceases  as  a function  of  applied  gap  voltage  for  various  indicated 
diameters  of  the  PTFE  insulator. 

a preamplifier,  A,  connected  to  an  impedance.  Z.  in  series  with  the  discharge  gap  as 
shown.  The  output  of  A was  fed  to  the  system  shown  in  figure  27 

The  system  was  used  here  to  measure  the  unconditional  pulse-height  and  pulse  time- 
separation  distributions,  po(qn)  and  po(-hfn)  respectively,  and  the  conditional  pulse- 
height  distribution  Pi  I Atn_i). 

It  is  found  for  the  configuration  shown  in  figure  33  that  there  is  a critical  point-to- 
plane  spacing  below  which  it  is  not  possible  to  observe  a recognizable  continuous 
Trichel-pulse  discharge  behavior.  Figure  34  shows  the  measured  critical  point-to- 
plane  gap  separations,  he,  at  which  the  Trichel-pulse  discharge  ceases  for  different 
indicated  values  of  the  diameter  D.  A pulsating  discharge  may  still  occur  for  d < dc, 
but  the  time  separation  between  pulses  becomes  very  long  (many  seconds).  It  is  seen 
that  the  larger  the  diameter  of  the  dielectric,  D,  and  the  lower  the  applied  voltage, 
14,  the  larger  will  be  the  critical  distance  d^  at  which  the  corona  extinguishes. 

Examples  of  measured  poiqn)  and  po(^^n)  at  different  indicated  gap  distances  d 
are  shown  in  figures  35  and  36  respectively  for  the  conditions  = 6.5kV,  D = 
5.0  cm  and  14  = 10.1  kV,  D = 7.5  cm.  For  convenience  in  plotting,  all  of  the 


Log  (Po(Atf,)/po(At 
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Figure  35.  Measured  pulse- height  and  corresponding  time-interval  distributions  at 
Va  = 6.5]lY,D  = 5.0cm,  for  the  various  indicated  gap  spacings.  The  distributions  have 
been  normalized  to  the  maxima. 
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Atn  iiLs) 

Figure  36.  Measured  pulse-height  and  corresponding  time-interval  distributions  at 
Va  = 10.1  kV,i?  = 7.5  cm,  for  the  various  indicated  gap  spacings.  The  distributions  have 
been  normalized  to  the  maxima. 
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Figure  37.  Mininnim  possible  pulse  time  separation,  Ate,  versus  gap  spacing  for 
Va  = 7.5  kV  and  the  different  indicated  PTFE  diameters. 

distributions  shown  in  these  figures  have  been  arbitrarily  normalized  to  the  maximum 
values.  The  profiles  of  the  distributions  shown  here  exhibit  features  similar  to  those 
previously  reported  and  explained  by  Van  Brunt  and  Kulkarni  [12].  For  example, 
the  distributions  po(Atn)  indicate  a previously  observed  [12, 72]  critical  minimum 
time  separation,  Ate,  between  successive  pulses.  The  rate  of  decrease  in  po(Atn) 
for  Atn  > At n, maxi  where  Atn,max  IS  the  maximum  in  po(Atn),  is  controlled  by  the 
rate  of  electron  emission  at  the  point  electrode  [12]  and  is  thus  strongly  dependent 
on  the  local  electric  field.  It  is  known  that  Ate,  which  is  controlled  by  the  ambipolar 
diffusion  of  ions  away  from  the  point,  also  decreases  exponentially  with  the  field  near 
the  cathode.  Examples  of  results  for  Ate  versus  d at  different  indicated  values  of  D 
are  shown  in  figure  37.  The  fact  that  both  Ate  and  the  width  of  go  through  a 

minimum  at  a particular  d > dc  indicates  that  the  electric  field  at  the  cathode  goes 
through  a maximum  at  this  separation.  This  maximum  in  the  field  does  not  occur 
when  the  dielectric  is  removed  as  illustrated  by  the  measured  data  for  po{qn)  and 
po{Atn)  shown  in  figure  38  that  were  obtained  under  conditions  where  no  dielectric 
was  present  on  the  planar  electrode.  The  corresponding  Ate  extracted  from  these  data 
are  plotted  in  figure  39  and  show  the  expected  monotonic  decrease  with  decreasing 
d which,  for  a fixed  Va,  is  consistent  with  an  increasing  field  strength  at  the  point 
electrode. 

The  profiles  of  the  pulse- height  distributions,  po{qn),  are  generally  consistent  with 
equation  (6)  and  thus  can  be  explained  in  terms  of  the  physical  processes  that 
determine  the  profile  of  the  corresponding  po{Atn)  distribution  coupled  with  the 
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Figure  38.  Measured  pulse-height  and  corresponding  time- interval  distributions  at 
Va  = 6.5  kV  for  various  indicated  gap  spacings  when  no  dielectric  is  on  the  anode  surface. 
The  distributions  have  been  normalized  to  the  maxima. 
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Figure  39.  Minimum  possible  pulse  time  separation,  versus  gap  spacing  for  the 

different  indicated  values  of  the  applied  voltage  Va  and  for  the  case  where  there  is  no 
dielectric  on  the  anode  surface. 
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known  dependence  of  on  Atn-i  as  previously  discussed  [12].  The  extent  to  which 
qn  depends  on  At„_i  is  determined  from  measurements  of  pi(^„lAt„_i). 

Examples  of  measured  conditional  distributions  pi(^„|At„_i)  at  different  indicated 
values  of  At„_i  and  d and  for  Va  = 10.1  kV  and  D = 7.5  cm  are  shown  in  figure  40 
together  with  the  corresponding  unconditional  distribution  poiqn)  indicated  by  the 
dashed  line.  A strong  positive  dependence  of  on  Atn-i  is  seen  from  these  data. 
This  behavior  is  consistent  with  data  shown  in  the  previous  section  and  can  again  be 
explained  by  the  effect  that  moving  negative-ion  space-charge  clouds  from  previous 
pulses  have  on  the  electric  field  at  the  cathode  when  the  next  discharge  pulse  develops. 

The  data  on  the  unconditional  distributions  such  as  those  shown  in  figures  35  and  36 
can  be  used  to  compute  the  expectation  values  {Qn)  and  (l/Atn)  corresponding 
respectively  to  mean  discharge-pulse  magnitude  and  repetition  rate.  Examples  of 
the  results  of  such  determinations  for  different  values  of  d,  D,  and  Va  are  shown 
in  figure  41.  Certain  general  trends  are  apparent  from  this  figure.  The  mean 
discharge-pulse  repetition  rate  generally  decreases  with  decreasing  Va , and  increasing 
D;  whereas  the  mean  discharge  magnitude  (Qn),  exhibits  the  opposite  behavior,  i.e., 
it  decreases  with  increasing  Va  and  decreasing  D.  Both  (1/At„)  and  (Qn)  niay  either 
increase  or  decrease  with  d depending  on  the  values  of  D and  d. 

Figures  42 and 43  show  plots  of  the  expectation  values  (^„(Atn_i))  versus  At„_i, 
where 

(7n(Atn— l))  ~ / (^n  I Atri— 1 ) . (l  — ) 

Jo 

These  results  which  were  derived  from  equation  (8)  using  data  on  the  conditional 
distributions  such  as  shown  in  figure  42,  indicate  the  strong  dependence  of  qn  on  Atn-i 
and  are  consistent  with  previous  results  [12]  except  at  short  distances  {d  < 0.6  cm) 
at  which  a change  in  slope  becomes  evident  as  At^-i  increases.  This  is  seen  most 
clearly  by  the  data  at  d = 0.4  cm  in  figure  43. 

Several  significant  conclusions  can  be  drawn  from  the  results  presented  here.  First, 
it  is  evident  from  the  measured  distributions  poiqn),  Po(At„),  and  pi(5n|At„_i)  that, 
except  for  point-to-dielectric  spacings  very  close  to  the  critical  value  dc  at  which 
Trichel-pulse  activity  ceases,  the  presence  of  the  dielectric  surface  has  little  or  no 
effect  in  modifying  the  general  stochastic  behavior  of  the  phenomenon.  For  sufficiently 
large  spacings,  the  observed  profiles  of  the  unconditional  pulse-amplitude  and  pulse- 
time separation  distributions  and  the  correlations  between  the  amplitudes  and  time 
separations  of  successive  pulses  are  consistent  with  those  reported  previously  [12]  for 
Trichel  pulses  that  occur  in  air  or  N2  — O2  mixtures  when  there  is  no  dielectric  on 
the  planar  electrode. 

It  can  be  argued  from  the  data  that  give  the  dependencies  of  Ate  and  (1/At„}  on  D 
shown  in  figures  37  and  41,  that  the  primary  effect  of  the  dielectric  on  Trichel-pulse 
development  at  any  spacing  d arises  from  a reduction  in  the  electric-field  magnitude  at 
the  point  electrode  caused  by  the  negative  charge  that  accumulates  on  the  dielectric 


54 


GASEOUS  DIELECTRICS  RESEARCH 


qn(pC) 

Figure  40.  Measured  conditional  pulse  amplitude  distributions  at  the 

indicated  values  for  Atn_i  and  gap  spacing,  and  the  corresponding  unconditional 
pulse- amplitude  distributions  po(?n)  shown  by  the  dashed  Hnes.  The  distribution  have 
been  normalized  to  the  maximum  values. 


1/(Atn)  {s-1)  <Qn)  (pC) 
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Figure  41.  Dependencies  of  the  mean  pulse  amplitude  (Qn)  a-nd  repetition  rate  l/(At„) 
on  gap  spacing  for  the  indicated  values  of  Va  and  D calculated  from  the  distributions  po{qn ) 
and  poi^^tn)  such  as  shown  in  figures  35 and 36. 
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Figure  42.  Tlie  expectation  value  (g„(At„_i))  versus  At„_i  at  = 10.1  kV,  D = 7.5  cm 
for  the  indicated  gap  spacings. 


Figure  43.  The  expectation  value  (g„(Atri-i))  versus  Atn-iforVa  = l0.1kV,D  = 7.5  cm 
and  d = 0.4  cm. 
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surface  during  the  discharge.  This  conclusion  is  supported  by  the  po(Atn)  data  in 
figures  35  and  38  for  d = 2.0  to  5.0  cm  that  were  obtained  respectively  with  and 
without  a dielectric  on  the  anode  surface.  The  observed  increases  in  both  Ate  and 
the  width  of  po(Atn)  which  occur  when  the  dielectric  is  introduced  are  consistent  with 
a reduction  of  the  electric-field  magnitude  equivalent  to  that  which  can  be  achieved 
by  reducing  the  gap  voltage.  As  previously  noted,  the  parameters  Afc,  {1/ At n),  and 
the  width  of  po{Atn)  are  sensitive  to  small  changes  in  electric-field  magnitude  at  the 
point  [12].  The  effect  of  increasing  the  diameter  of  the  dielectric  surface  is  also  similar 
to  that  of  decreasing  the  gap  voltage.  The  appearances  of  corresponding  minima  in 
Ate  and  the  widths  of  the  po{Atn)  distributions  at  a particular  gap  spacing  indicate 
that  the  field  at  the  point  electrode  must  attain  a maximum  at  that  spacing  for  a 
fixed  Va  when  the  dielectric  is  present. 

It  can  be  seen  from  figures  42  and  43  that  the  presence  of  the  dielectric  begins  to 
influence  the  dependence  of  on  Atn-i  significantly  for  small  point-to-dielectric  gap 
spacings  close  to  the  critical  value  dc.  From  an  examination  of  the  data  shown  in 
figure  43,  it  can  be  argued  that  for  A^n-i  < Af^,  the  dependence  of  on  Atn-i 
is  due  predominantly  to  the  previously  mentioned  negative-ion  space-charge  effect 
which  is  a prominent  characteristic  in  the  stochastic  behavior  of  Trichel  pulses  in  the 
absence  of  the  dielectric  barrier  [12].  The  slope  of  the  (g„(Af„_i))  versus  Af„_i  curve, 
however,  is  seen  to  decrease  significantly  at  At^.  It  is  speculated  that  at  Atn-i  = Ati, 
the  negative  ions  clear  the  gap  and  are  presumably  deposited  on  the  dielectric.  Thus, 
for  Atn-i  > Ai^,  the  dependence  of  on  A^„_i  could  be  governed  primarily  by  the 
rate  of  charge  dissipation  or  migration  on  the  dielectric  surface.  It  should  be  noted 
that  the  additional  surface  charge  deposited  by  an  earlier  discharge  pulse  which  is 
present  when  the  next  pulse  initiates  will  act  to  reduce  the  field  at  the  point  electrode 
when  this  pulse  develops  and  will  thereby  tend  to  inhibit  its  growth.  Therefore,  a 
dissipating  residual  surface  charge  from  a previous  pulse  acts  like  the  moving  residual 
negative-ion  space  charge  to  cause  to  have  a positive  dependence  on  Atn_i. 

Measurements  performed  with  a vibrating  plate  type  electric-field  probe  [73,74] 
showed  that  after  the  corona  discharge  is  turned  off,  the  PTFE  retains  a quasi- 
permanent negative  surface  charge  that  decays  slowly  with  a time  constant 
(>  lO'*  seconds)  which  is  much  greater  than  the  mean  time  between  Trichel  pulses. 
This  observation  suggests  that  if  the  dependence  of  on  A^„_i,  for  A^„_i  < At(  in 
figure  43  is  indeed  due  to  an  effect  of  surface- charge  dissipation,  it  must  correspond  to 
a component  of  the  charge  density  that  is  much  more  mobile  than  that  detected  with 
the  field  probe,  i.e.,  it  must  decay  with  a time  constant  comparable  in  magnitude  to 
the  mean  time  between  discharge  pulses.  This  is  consistent  with  recent  observations 
of  Jonscher  and  Owede  [75]  which  indicate  that  the  charge  deposited  on  insulating 
surfaces  can  become  increttsingly  more  mobile  as  the  charge  density  is  increased. 

The  conclusions  drawn  here  about  the  influence  of  a PTFE  dielectric  surface  on 
Trichel-pulse  discharge  can  be  expected  to  apply  in  general  to  other  types  of 
insulating  materials  to  the  extent  that  these  materials  can  acquire  permanent  or 
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quasi-permanent  charge  from  the  corona.  A preliminary  survey  made  using  other 
materials  suggested  that  this  is  indeed  the  case. 


2.3.5  Limitations  and  Extensions  of  the  SAPP  Technique 

The  system,  as  shown  in  figure  27,  has  certain  limitations.  First,  it  is  designed  to 
operate  with  dc  voltages.  Second,  as  discussed  in  ref.  [60],  it  cannot  give  reliable 
measurements  of  certain  conditional  distributions  if  the  time  separations  between 
events  become  shorter  than  either  the  TAG  reset  time  or  the  MCA  analog- to- digital 
conversion  time.  Third,  it  cannot  be  applied  to  nonstationary  phenomena  that  drift 
on  a time  scale  comparable  to  the  time  required  to  acquire  statistically  significant 
data.  Fourth,  it  can  only  measure  one  distribution  at  any  given  time.  Fifth, 
independent  of  the  problems  associated  with  timing  in  either  the  TAG  or  MCA,  it 
does  not  employ  broad-band  detection  methods  and  therefore  cannot  resolve  pulses 
separated  by  times  less  than  the  width  of  the  impulse  response.  Because  of  this 
limitation,  it  could  not,  for  example,  be  applied  to  investigate  P-D  pulse  bursts 
which  as  shown  below  occur  in  liquids  and  are  known  to  occur  in  the  positive  corona 
of  some  highly  electronegative  gases  like  SFs  [55]. 

The  restriction  to  dc  voltages  can  be  overcome  by  incorporating  a phase  marker  such 
as  a zero-crossing  sensor  to  identify  the  phase  associated  with  any  given  pulse  time 
separation.  The  system  has,  in  fact,  been  recently  modified  to  allow  measurement 
of  conditional  distributions  restricted  to  specific  phase  intervals.  For  example  it  now 
allows  measurement  of  the  phase  conditioned  distributions,  p\{qi((pi)  \ where 

qi  is  the  amplitude  of  any  P-D  pulse  that  occurs  with  phase  (pi  contained  in  the 
specified  interval  /\(pj.  This  type  of  distribution  has  been  considered  recently  by 
others  [58, 59].  However,  in  addition  to  this  distribution,  our  system  is  capable  of 
measuring  many  others  in  which  phase  is  specified.  These  include  the  distributions 
Pi{qi(i\(t>2)  I Q{/A(f>i))  and  pi(5i(A(^2  ) I where  q\  is  the  amplitude  of  the 

first  pulse  to  occur  in  either  the  positive  or  negative  half  cycle  specified  respectively 
by  the  phase  intervals  A(/>2  or  A<2!>7  and 

Q((\<t>t)  = <pf  € Aaif,  (13) 

is  the  total  charge  associated  with  all  P-D  pulses  occurring  within  phase  interval 
of  the  preceding  half  cycle.  Preliminary  data  obtained  on  these  distributions 
for  dielectric-barrier  discharges  at  60-Hz  voltages  show  that  the  charge  deposited  on 
a dielectric  surface  by  partial  discharges  during  one  half  cycle  significantly  affect  the 
development  of  partial  discharges  during  the  subsequent  half  cycle. 

Limitations  on  minimum  time  separations  can  be  overcome  in  principle  by  using 
prerecorded  data  that  can  be  obtained  with  fast-transient  digitizers.  However,  if 
prerecorded  data  are  to  be  used,  it  may  be  preferable  to  use  computer  software  to 
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determine  the  conditional  distributions,  rather  than  the  hardware  presented  here.  It 
should  be  emphasized  that  an  important  advantage  of  the  method  described  here 
is  that  the  conditional  distributions  are  measured  directly  in  real  time.  This  is  an 
advantage  because  it  overcomes  the  problem  of  having  to  store  enormous  quantities 
of  data  required  to  determine  conditional  distributions  by  the  software  approach.  For 
example,  some  of  the  second-order  conditional  distributions  shown  in  figures  28  and  29 
required  about  10  minutes  of  data  acquisition  time  in  order  to  obtain  acceptable 
statistics.  This  means  that  data  on  all  but  a very  small  fraction  of  the  discharge  events 
that  occurred  in  this  10  minute  period  were  discarded  by  the  filtering  procedure. 
Using  the  software  approach,  all  the  pulses  that  occurred  during  that  time  would 
have  to  be  stored,  which  for  the  data  shown  means  that  more  than  10^  pairs  of 
numbers  would  have  to  be  acquired  just  to  determine  one  distribution.  Also,  because 
the  present  system  works  in  real  time,  the  distributions  can  be  seen  immediately  as 
the  data  are  accumulated. 

If  the  P-D  phenomenon  is,  or  can  become,  nonstationary,  then  an  independent  method 
may  be  needed  to  determine  the  extent  to  which  it  is  nonstationary.  A periodic  check 
of  the  unconditional  pulse- amplitude  and/or  time-separation  distribution  will,  in  most 
cases,  provide  an  indication  of  deviations  from  stationary  behavior. 

In  the  case  of  multiple  P-D  discharge  sites,  an  ambiguity  can  occur  in  the 
interpretation  of  data  from  measurement  of  conditional  distributions,  particularly 
if  the  discharges  from  the  different  sites  occur  with  nearly  equal  intensity.  The 
measurement  scheme  discussed  here  works  best  if  there  is  only  one  site  at  which 
P-D  discharges  are  being  generated  with  an  intensity  that  is  significantly  above  the 
background  P-D  level.  The  existence  of  pronounced  correlations  among  successive 
P-D  pulses  would,  in  fact,  indicate  that  the  discharges  are  occurring  at  one.  or  a 
limited  number  of  sites.  The  determination  of  this  fact  might  be  important  in  some 
diagnostic  applications. 

Finally,  it  should  be  pointed  out  that  the  present  system  could  be  modified  to  allow  a 
more  efficient  use  of  the  available  data.  One  could,  for  example,  use  multiple  MC.A.’s 
in  parallel  to  measure  simultaneously  more  than  one  distribution. 
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Task  03 
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3.1  Introduction 

Electric  breakdown  of  liquid  dielectrics  is  preceded  by  the  growth  of  a streamer 
in  the  liquid.  Recent  studies  of  these  phenomena,  see  for  example  Hebner[76]  for 
a review,  reveal  that  as  a streamer  propagates  it  follows  a well  defined  pattern: 
Each  stage  or  mode  of  growth  may  be  distinguished  by  its  structure,  speed  of 
propagation,  and  relative  extent.  Considerable  effort  has  been  devoted  to  describing 
the  conditions  that  influence  streamer  propagation  and  initiation  [77].  However, 
despite  this  activity,  no  clear  consensus  has  been  reached  as  to  the  mechanism 
for  the  initiation  of  the  discharge.  Proposals  put  forward  have  included:  electron 
avalanche  [78],  cavitation  [79],  and  electrostatic  forces  acting  on  the  dielectric  [80]. 

The  mechanisms  responsible  for  the  initial  growth  of  negative  streamers  appear  to 
be  distinct  from  those  during  its  latter  stages.  Watson  and  Chadband[81]  show 
initial  growth  of  negative  streamers  to  be  consistent  with  the  propagation  of  a cavity 
within  the  liquid.  Similar  results  are  reported  for  partial  discharges  (PD’s),  where 
photographic  and  optical  studies  [78,  82,  83]  have  established  a clear  correspondence 
between  the  growth  of  a cavity  in  the  fiquid  and  bursts  of  discrete  current  pulses 
emitted  from  a cathode.  The  results  presented  in  this  report  were  obtained  during 
experimental  investigations  of  the  initiation  and  early  growth  of  negative  streamers  in 
hexanes  [84,  85].  In  section  3.2,  the  initial  growth  of  negative  streamers  generated  by 
an  impulse  voltage  is  examined  at  pressures  ranging  from  0.1  to  1.1  MPa.  The  pressure 
dependence  of  the  streamer  initiation  and  inception  voltages  are  shown.  In  section  3.3, 
partial  discharges  generated  by  the  application  of  dc  voltages  are  examined.  A non- 
uniform  field  geometry  is  employed  and  the  growth  of  cavities  associated  with  partial 
discharges  at  a point  cathode  are  photographed  at  high  magnification.  A simultaneous 
record  of  the  discharge  current  is  also  obtained.  This  study  employs  an  image- 
preserving optical  delay  [82,  82],  and  the  use  of  this  device  allows  photography  of 
the  discharge  at  its  inception.  Such  data  allow  detailed  description  of  the  temporal 
and  spatial  development  of  the  cavity  and  provide  a basis  for  the  evaluation  of  models 
for  the  initiation  of  negative  streamers. 
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3.2  Streamer  Initiation  in  Hexanes 

For  hydrocarbon-based  dielectric  liquids,  four  modes  of  streamer  growth  are 
suggested  [76, 87]  for  negative  streamers;  these  are  shown  schematically  in  figure  44. 
The  1st  mode  appears  to  be  a single,  relatively  linear  channel  emanating  from  the 
cathode  surface  and  connecting  directly  to  a much  more  bushy  structure,  mode  2. 
During  modes  1 and  2,  the  streamer  propagates  at  subsonic  speeds.  Although  the 
3rd  mode  originates  at  the  2nd  mode,  its  structure  differs  from  that  of  the  2nd  mode, 
and  its  speed  of  propagation  is  significantly  higher,  approaching  the  sonic  velocity. 
The  transition  from  the  2nd  to  the  3rd  mode  is  not  continuous;  streamer  growth 
often  stalls  after  the  formation  of  the  2nd  mode  [88].  The  4th  mode  directly  precedes 
breakdown  and  it  appears  to  be  a single  streamer  that  propagates  at  supersonic 
speeds. 

Experimental  conditions  influence  streamer  propagation;  for  example,  the  relative 
extent  of  the  2nd  mode  is  sensitive  to  both  the  ambient  pressure  [89-91]  and  the  purity 
of  the  fluid  [92].  Indeed,  it  has  been  suggested  that  at  sufficiently  high  pressures,  the 
1st  and  2nd  modes  may  not  exist  at  all  [89,  90,  93].  Further  systematic  study  of  these 
effects  is  of  value;  and  the  present  work  was  undertaken  to  examine  the  inception 
and  early  growth  of  negative  streamers.  For  this  purpose,  highly  magnified  frame 
photographs  of  streamer  growth  at  a needle  tip  as  well  as  photographs  of  the  full  gap 
are  obtained  at  pressures  ranging  from  0.1  - 1.1  MPa. 

3.2.1  Apparatus  and  Procedure 

Streamer  growth  is  photographed  using  conventional  shadowgraphic  techniques  and 
an  image-converter  camera  [87, 89].  The  optical  resolution  is  improved  over  that 
obtained  in  previous  studies  by  mounting  a microscope  objective  within  the  test 
cell.  This  configuration  provides  an  optical  resolution  of  approximately  1 ^m  at 
100  X magnification. 

The  test  gap  consists  of  a steel  gramophone  needle  and  a stainless  steel  rod;  the 
electrodes  are  enclosed  within  a brass  cell  and  are  immersed  in  the  test  fluid.  The 
curvature  of  the  needle  is  approximately  parabolic,  having  a radius  of  curvature  at 
its  tip  of  approximately  40 //m.  The  rod  electrode  is  6.4  mm  in  diameter  and  has  a 
hemispherical  tip.  The  needle  and  rod  are  separated  by  a distance  of  5 mm. 

A voltage  impulse,  having  negative  polarity,  is  applied  directly  to  the  needle  and 
the  rod  electrode  is  connected  to  ground.  The  voltage  waveform  is  approximately 
trapezoidal:  the  rise  and  fall  times  are  3/zs  (10-90%)  and  the  voltage  varies  by  less 
than  ±10%  within  a 4//s  interval  centered  about  the  peak  voltage.  .A.  precision 
resistive  voltage  divider  and  digital  oscilloscope  are  used  to  record  the  waveform. 
The  camera  monitor  pulse  is  also  recorded  and,  together  with  the  voltage  waveform, 
these  data  determine  the  voltage  at  each  frame  within  a photograph.  Two  measures 
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of  the  effects  of  pressure  are  determined  by  these  methods:  the  average  streamer 
inception  voltage,  Vine,  and  initiation  voltage,  Vinit-  For  purposes  of  this  study,  Vinit 
is  the  lowest  voltage  for  which  one  streamer  is  observed  in  eight  impulses,  and  is  a 
measure  of  the  minimum  voltage  required  to  produce  a prebreakdown  streamer.  For 
voltages  greater  than  or  equal  to  the  inception  voltage,  a second  useful  quantity  may 
be  determined:  the  average  voltage  at  the  inception  of  the  streamer  Vine-  Values 
of  Vinit  and  Vine  are  measured  using  high  magnification;  furthermore,  to  avoid  rapid 
deterioration  of  the  needle  tip  and  of  the  liquid,  they  are  obtained  without  breakdown. 

Ultra-pure  grade  hexanes  are  used  both  for  this  study  and  for  the  work  discussed  in 
section  3.3.  The  liquid  is  used  as  received  without  degassing  or  purification  other  than 
by  filtration.  While  at  atmospheric  pressure,  the  liquid  is  continuously  circulated 
and  filtered  through  a 2 /im  pore  size  paper  filter.  The  pressure  within  the  cell 
is  increased  by  compressing  the  liquid,  which  is  accomplished  by  the  use  of  small 
hydraulic  cylinder  connected  to  the  cell.  The  application  of  pressure  did  not  allow 
continuous  filtration;  however,  pressure  was  released  and  the  liquid  circulated  and 
filtered  every  8-10  impulses. 

3.2.2  Results  and  Discussion 

Photographic  Observation  of  Streamer  Behavior. 

Photographs  of  the  full  gap  obtained  at  pressures  of  0.1,  0.7,  and  1.1  MPa  are  shown 
in  figure  45.  The  effects  of  pressure  are  made  evident  in  these  photographs  by  a 
dramatic  reduction  in  the  extent  of  the  2nd  mode  and  by  changes  in  the  structure  of 
the  3rd  mode:  it  may  be  noted  that  the  latter  is  less  dense  and  has  fewer  branches  at 
high  pressure.  The  1st  and  2nd  modes  may  be  readily  identified  in  the  photographs 
obtained  at  atmospheric  pressure  and  remain  visible  up  to  0.7  MPa,  however  they 
are  not  resolved  in  the  low-magnification  photographs  for  pressures  of  0.9  MPa  and 
above. 

For  comparison,  high-magnification  photographs  obtained  at  pressures  of  0.1,  0.5, 
and  1.1  MPa  are  shown  in  figure  46.  In  contrast  to  the  photos  taken  at  low 
magnification,  the  characteristic  stem  and  bush  structure  is  clearly  evident  at  all 
pressures.  Furthermore,  the  images  are  sufficiently  well-resolved  to  indicate  the 
dynamical  behavior  of  the  streamer:  upon  close  examination  of  these  photographs, 
it  is  evident  that  the  1st  mode  expands  and  contracts  repeatedly  as  the  streamer 
propagates.  It  was  further  noted,  that  the  frequency  and  amplitude  of  oscillation  of 
the  1st  mode  increase  with  pressure.  Indeed,  at  the  highest  pressure  examined,  the 
1st  mode  completely  disappears  and  reappears.  The  maximum  diameter  and  length 
of  the  1st  cathode  mode  are,  however,  relatively  insensitive  to  pressure  and  are  on 
the  order  4 pm  and  S pm  respectively.  The  2nd  mode  did  not  exhibit  such  dramatic 
oscillations. 
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The  oscillations  in  the  diameter  of  the  1st  cathode  mode  may  be  consistent  with 
intermittent  light  emission  and  charge  injection  observed  during  the  propagation  of 
negative  streamers  [94,  95].  Studies  of  partial  discharge  phenomena  at  a point  cathode 
in  the  same  material  [85]  show  current  to  be  highly  correlated  with  the  presence  of  a 
bubble- like  structure  at  the  electrode.  It  seems  likely  that  significant  charge  injection 
can  only  occur  when  the  streamer  is  in  contact  with  the  electrode;  and  thus  that  the 
prebreakdown  streamer  provides  a preferred  path  for  current  flow.  Indeed,  studies  of 
full  breakdown  show  the  conducting  channel  to  lie  along  the  path  of  a prebreakdown 
streamer. 

Streamer  Initiation  and  Inception  Voltages. 

The  methods  described  allow  the  effects  of  changes  in  the  ambient  pressure  and  in 
the  peak  voltage  on  the  streamer  inception  voltage  to  be  examined.  Values  of  Vine 
are  estimated  from  eight  trial  pulses  and  are  measured  using  a 217  ns  frame  interval. 
The  results  of  these  measurements  for  pressures  of  0.1, 0.5,  and  0.9  MPa  are  shown  in 
figure  47,  where  Vine  are  plotted  on  representative  traces  of  the  voltage  waveform. 
Clear  trends  in  both  the  formation  time  and  inception  voltage  are  apparent.  At 
0.1  MPa,  Vine  has  a strong  dependence  on  the  voltage  waveform:  the  time  required 
for  streamer  inception  decreases  and  the  voltage  at  inception  increases  with  applied 
voltage.  Although  the  formation  time  decreases  with  faster  rising  impulses,  it  may 
be  noted  that  the  inception  voltage  is  relatively  constant  for  the  data  obtained  at 
0.9  MPa.  Indeed  Vine  appears  to  approach  a value  that  is  relatively  independent  of 
pressure  as  the  peak  voltage  is  increased.  These  observations  lend  support  to  the 
results  of  Kao  and  McMath[96],  who  show  for  both  transformer  oil  and  n-hexane 
that  increased  pressure  has  little  or  no  influence  on  breakdown  voltage  for  fast  rising 
voltages.  In  their  study,  the  effect  of  increased  pressure  was  found  to  be  negligible 
when  the  rate  of  rise  of  the  applied  voltage  exceeded  4MV/cm/is.  For  the  results 
presented  here,  the  rate  of  rise  increases  with  applied  voltage.  Numerical  methods  [97] 
used  to  estimate  the  electric  field  at  the  needle  tip  yield  5.5  x 10^  V/cm  at  a 30  kV 
applied  potential.  This  value  together  with  the  greatest  rate  of  voltage  rise  shown  in 
the  figure  yield  6.6MV/cm/is  which  compares  favorably  with  the  value  given  by  Kao 
and  McMath. 

The  effect  of  pressure  on  the  streamer  initiation  voltage  is  shown  in  figure  48.  Vi^it  is 
seen  to  increase  with  pressure  from  0.1  to  0.5  MPa  and  to  be  relatively  independent 
of  pressure  above  0.5  MPa.  Peak  voltages  of  —33  and  — 46kV  were  used  at  0.1  and 
0.3  MPa,  respectively,  and  — 53kV  for  pressures  of  0.5  MPa  and  above.  In  each  case, 
streamer  initiation  occurred  near  the  peak  voltage. 
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1 St  Cathode  Mode 


Figure  44.  Cathode  streamer  development.  The  four  modes  of  streamer  growth  are  shown, 
(Hebner  [76]). 
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(c)  1.1  MPa 


Figure  45.  Low  magnification  photographs  of  prebreakdown  streamers.  The  needle 
electrode  is  shown  on  the  right  in  each  of  these  photographs.  Modes  1 and  2 may  be 
readily  identified  in  the  photographs  obtained  at  0.1  MPa.  The  frame  interval  is  512  ns  and 
the  frame  sequence  is  shown  in  (a).  An  impulse  having  a peak  voltage  of  —93  kV  is  applied 
to  the  needle. 
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Figure  46.  High- magnification  photographs  of  prebreakdown  streamers.  The  pressure  and 
peak  voltages  axe  indicated  beside  each  photograph  and  the  frame  interval  is  217  ns. 
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Figure  47.  Pressure  and  waveform  effects  on  streamer  inception.  Tiie  average  streamer 
inception  voltages  Vine,  a-re  plotted  on  representative  traces  of  the  voltage  waveform.  Values 
obtained  at  pressures  of  Q,  0.1;  A,  0.5;  and  □,  0.9  MPa  are  shown. 
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Figure  49.  Partial  discharge  record  obtained  at  17.5  kV.  The  frame  interval  and  exposure 
are  51  and  10  ns  respectively.  The  current  waveform  and  the  extent  of  the  cavity  are  plotted 
in  the  upper  panel. 
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3.3  DC  Voltage  Conditions 

For  dc  voltages  near  the  threshold  for  P-D  activity,  P-D’s  occur  infrequently  and 
are  separated  by  time  intervals  that  are  widely  distributed.  Furthermore,  significant 
current  may  only  persist  for  times  of  the  order  of  microseconds.  These  characteristics 
present  a technical  challenge:  photographing  a short  lived  event  that  may  occur 
essentially  at  random.  The  methods  employed  allow  two  independent  records  of 
the  P-D  to  be  obtained  simultaneously:  measurement  of  P-D  current  and  frame 
photography  of  the  cavity. 

3.3.1  Experimental  Apparatus 

The  test  gap  consists  of  a steel  needle  and  stainless  steel  rod  separated  by  3.2  mm. 
The  electrodes  are  enclosed  within  a brass  cell  and  are  immersed  in  the  test  liquid. 
The  liquid  is  used  as  received  without  degassing  or  further  purification.  The  rod 
electrode  is  6.4mm  in  diameter  and  has  a hemispherical  tip,  and  is  connected  to  a 
high-voltage  dc  power  supply.  The  needle  electrode  is  1.0  mm  in  diameter  and  has  a 
conical  point  with  an  apex  angle  of  30°  and  a tip  radius  of  1 pm. 

A continuous- wave  argon  laser  is  used  to  illuminate  the  needle  tip,  and  frame 
photographs  of  the  discharge  are  obtained  by  shadowgraphic  methods  and  the  use 
of  an  image- converter  camera.  The  total  magnification  is  200  x and  the  optical 
resolution  is  sufficient  to  fully  resolve  objects  less  than  2pm  in  diameter.  The  use  of 
an  image-preserving  optical  delay  allows  the  framing  sequence  to  begin  much  closer 
to,  and  for  some  frame  intervals  before,  the  on  set  of  the  partial  discharge. 

A broadband  transimpedance  amplifier  [98,  99]  is  used  to  measure  the  P-D  current. 
The  amplifier  is  designed  to  provided  a low-noise  background  against  which  low- 
level  P-D  currents  may  be  detected.  The  amplifier  has  three  stages:  a low-noise, 
broad-bandwidth,  current-sensitive  preamplifier  and  low-gain  and  high-gain  buffer 
amplifiers.  The  combined  system  bandwidth  is  estimated  to  be  35  MHz  and  the 
equivalent  noise  is  30nA(rms).  The  charge  sensitivity  of  the  system  is  estimated  to 
be  less  than  0.7  fC.  The  amplifier  is  connected  directly  to  the  needle  electrode  and  the 
current  waveform  is  recorded  digitally.  A digital  oscilloscope  having  a sample  rate  of 
10®  samples/s  is  used  to  obtain  the  results  shown  here.  Both  the  digital  oscilloscope 
and  the  framing  camera  derive  their  trigger  signals  from  the  P-D  current. 

3.3.2  Partial  Discharge  Records 

An  example  of  a P-D  record,  obtained  at  an  applied  potential  of  17.5  kV,  is  shown  in 
figure  49.  Note  that  the  first  few  camera  frames  precede  the  initiation  of  the  current 
waveform.  The  data  at  hand  allow  an  estimate  of  the  time  interval  between  the 
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initial  current  pulse  and  the  onset  of  cavity  growth.  As  shown,  the  current  waveform 
begins  51  ± 5 ns  before  the  midpoint  of  the  sixth  frame,  the  first  frame  in  which 
clear  evidence  of  the  cavity  appears.  The  estimated  uncertainty  in  the  position  of 
frame  six  is  ±8  ns;  combining  these  uncertainties  in  quadrature  yields  51  ± 9 ns  and 
thus  provides  an  upper  bound  of  60  ns.  The  extent  of  the  cavity  in  the  sixth  frame, 
however,  suggests  that  the  time  interval  is  likely  to  be  significantly  less.  In  the  upper 
panel,  the  length  of  the  cavity  is  plotted  against  the  frame  number  and  fit  to  a linear 
function.  Such  a fit  provides  an  estimate  for  the  origin  of  the  cavity  that  is  nearly 
coincident  with  the  initiation  of  the  current  waveform.  However,  in  the  absence  of  a 
more  accurate  determination  of  the  instantaneous  rate  of  expansion  of  the  cavity  at 
its  inception,  we  are  reluctant  to  apply  this  correction  and  prefer  to  report  the  upper 
bound. 

The  correlation  between  the  onset  for  cavitation  and  PD  current  noted  here  lends 
support  to  recent  experimental  results  obtained  in  cyclohexane  [78].  Results  obtained 
in  nitrobenzene  suggest  that  significant  charge  may  be  injected  into  the  liquid  prior  to 
streamer  initiation  [100,  101],  and,  that  the  onset  for  charge  injection  precedes  that 
for  cavitation  by  as  much  as  4/zs[100].  The  presence  of  charge  in  the  fluid  would 
result  in  a force  on  the  fluid  of  the  proper  sign  to  produce  the  observed  cavitation. 

The  linear  fit  also  provides  an  estimate  of  the  average  rate  of  expansion  of  the  cavity. 
Furthermore,  since  the  rate  of  expansion  depends  on  the  forces  acting  on  the  liquid,  an 
estimate  of  the  effective  pressure  within  the  cavity.  Pin,  may  be  obtained.  Assuming 
inviscid  incompressible  flow,  and  approximating  the  cavity  by  an  expanding  sphere 
of  radius  r,  the  instantaneous  work  done  in  expanding  the  cavity  is 


25  3 2 

Pin  dv  = Pamb  dv  -{ dv  + -pr  dv. 

r 2 


The  terms  include:  the  work  done  against  ambient  pressure,  Pamb?  aJid  surface  tension, 
5;  and  the  kinetic  energy  imparted  to  the  surrounding  liquid.  The  average  rate  of 
expansion  of  the  cavity,  f,  is  approximately  20  m/s  and  the  density,  /?,  and  surface 
tension  of  hexanes  are  approximately  0.66  g/cm^  and  1.8  x 10”^N/m  respectively. 
For  the  conditions  described,  the  kinetic  energy  is  the  largest  term  and  P\n  is 
approximately  0.4  MPa. 

The  required  pressure  may  be  provided  by  local  heating  of  the  liquid,  as  may  occur 
during  electron  avalanche,  and  by  electrostatic  forces.  The  form  of  the  cavity, 
however,  suggests  that  the  forces  acting  on  the  liquid  are  highly  directional,  and  thus 
that  the  gas  dynamic  pressure  within  the  cavity  is  of  limited  importance.  We  believe 
these  observations  lend  support  to  the  results  reported  by  Watson  and  Chadband  [81] 
and  further  suggest  that,  from  its  inception,  the  growth  of  the  cavity  may  be  described 
by  the  electrostatic  forces  acting  on  the  liquid. 

The  current  waveform  is  relevant  to  describing  the  conditions  at  the  initiation  of  the 
P-D.  In  over  90%  of  the  current  waveforms,  clear  evidence  of  a continuous  current 
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is  detected  at  the  initiation  of  the  waveform.  Furthermore,  whereas  the  later  current 
pulses  tend  to  grow  monotonically  in  amplitude,  the  first  two  pulses  are  frequently  of 
comparable  amplitude.  Indeed,  P-D’s  are  observed  that  consist  of  only  a double- 
pulse structure.  These  observations  suggest  that  the  mechanisms  producing  the 
measured  current  at  the  origin  of  the  cavity  growth  differ  from  those  at  later  times. 
Two  mechanisms  that  give  rise  to  continuous  currents  are  suggested:  the  motion  of 
charge  trapped  in  the  liquid  or  at  the  interface  between  the  liquid  and  the  metal 
electrode  would  produce  a current  in  the  external  circuit  [102],  also  the  reduction  in 
the  capacitance  of  the  cell  due  to  a reduction  in  the  dielectric  constant  of  the  material 
near  the  needle  tip  would  produce  a continuous  current. 
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3.4  Conclusions 

The  effects  of  increased  ambient  pressure  on  the  initiation  and  growth  of  negative 
streamers  in  hexanes  are  observed  using  a 100  x high-magnification  optical  system 
with  1 /zm  resolution.  The  existence  of  the  1st  cathode  streamer  mode  is  clearly  shown 
at  pressures  ranging  from  0.1  to  1.1  MPa.  The  1st  cathode  mode  is  approximately 
4 fim  in  diameter  and  8 /zm  in  length  at  the  transition  to  the  2nd  mode.  Dynamical 
behavior  of  the  1st  cathode  mode  is  shown;  the  diameter  of  the  stem  appears  to 
oscillate  as  the  streamer  grows  and,  at  high  pressure,  the  stem  completely  disappears 
and  reappears  although  the  2nd  mode  continues  to  exist.  The  behavior  of  the 
1st  cathode  mode  may  be  associated  with  intermittent  charge  injection  and  the 
stepped  propagation  of  negative  streamers  noted  in  previous  studies.  The  3rd  mode 
becomes  thinner  and  less  branched  at  higher  pressure.  The  streamer-initiation  voltage 
increases  with  pressure,  but  the  effect  of  pressure  on  the  streamer-inception  voltage 
decreases  at  higher  applied  voltage. 

Highly  resolved  photographs  of  the  initiation  of  P-D’s  in  liquid  hexanes  are  obtained 
for  dc  voltages.  The  cavity  growth  at  a point  cathode  is  nonisotropic,  which  suggests 
that  electrostatic  forces  are  of  primary  importance  in  driving  its  expansion.  The  onset 
for  P-D  current  precedes  or  is  simultaneous  with  the  initial  growth  of  the  cavity.  An 
upper  bound  of  60  ns  between  the  first  current  pulse  and  the  growth  of  the  cavity  is 
estimated. 
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4.1  Introduction 

The  purpose  of  an  impulse  voltage  measuring  system  is  to  provide  a means  of 
reducing  the  high  voltage  signal  to  levels  which  are  compatible  with  data  recording 
equipment  and  this  requires  that  the  scale  factor  of  the  system  be  constant 
throughout  the  frequency  range  of  interest.  Most  impulse  measurement  systems 
used  by  industrial  laboratories  invariably  introduce  some  distortion  of  the  input 
signal  due  primarily  to  inadequate  bandwidth,  voltage  coefficient,  and  some  other 
factors.  This  distortion  may  be  either  negligible  or  totally  unacceptable,  depending 
on  the  allowable  error  associated  with  the  particular  measurement  requirement. 
According  to  lEC  and  IEEE  standards  [103, 104],  a system  which  is  used  to  measure 
standard  lightning  impulses  should  have  errors  of  less  than  3 % and  10  % for  peak 
voltage  and  front  time,  respectively.  The  standards  also  recommend  that  certain 
aspects  of  the  dynamic  behavior  of  the  measurement  system  should  be  evaluated  by 
using  parameters  obtained  from  the  step  response.  There  is,  however,  no  reliable 
relationship  between  the  errors  of  a measurement  system  and  its  step  response 
parameters.  In  fact,  it  has  been  shown  that  some  measurement  systems  exhibit 
step  response  parameters  that  fall  within  the  limits  given  by  lEC  and  IEEE  but 
that  actually  produce  unacceptably  large  errors  in  the  measurement  of  high  voltage 
impulses  [105, 106]. 

A more  reliable  and  simplified  method  has  been  recommended  in  the  new  draft  of 
IEC-60  that  is  based  on  simultaneous  measurements  of  a high  voltage  impulse  by 
an  independent  reference  system  and  the  system  under  test  [107].  As  a first  step, 
international  comparative  measurements  were  made  in  four  National  Laboratories 
and  the  relative  differences  among  them  were  reported  [109].  Investigations  have 
also  been  made  of  the  interactions  between  two  systems  configured  for  simultaneous 
mecisurements,  and  of  methods  for  minimizing  these  interactions  [106]. 

The  comparative  method  establishes  the  uncertainty  in  measurements  of  a particular 
type  of  impulse  for  a measurement  system  by  comparison  against  a reference 
mecLSurement  system  having  overall  uncertainties  for  peak  voltage  of  less  than  1 % 
over  its  range  of  use  [107].  The  key  to  the  comparative  measurement  technique  is  that 
the  reference  system  has  to  be  built  and  carefully  characterized.  Without  a stable 
and  precise  reference  system,  it  is  usually  impossible  to  use  the  comparative  method 
to  estimate  the  uncertainty  of  another  system.  There  are,  however,  techniques  that 
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utilize  a pair  of  complementary  measurement  systems  to  detect  measurement  errors  in 
each  system,  but  the  technique  requires  that  the  measurement  systems  have  different 
error  sources  [108]. 

As  a part  of  a series  of  investigations  into  high  voltage  impulse  measurements,  a 
resistive  impulse  measurement  system  was  specially  designed,  built,  and  calibrated 
at  NIST.  It  is  cinticipated  that  this  system  along  with  some  Kerr  electro-optical 
devices  will  be  used  as  the  NIST  reference  systems.  This  report  documents  the 
design,  construction,  and  characterization  of  this  fast  impulse  voltage  measurement 
system,  which  is  rated  for  peak  voltages  up  to  300  kV. 


4.2  Resistive  Voltage  Divider  Design 

4.2.1  General  Considerations 

The  resistive  voltage  divider,  designated  as  NIST4,  serves  as  the  primary  component 
of  the  reference  measurement  system.  It  was  designed  to  accurately  measure 
full  lightning  impulses  up  to  300  kV.  The  voltage  range  was  selected  for  several 
reasons.  First,  it  was  similar  to  devices  used  as  reference  systems  in  the  national 
laboratories  of  other  countries  so  that  direct  comparison  between  the  devices  would 
be  reasonable.  Second,  compared  with  a divider  rated  for  megavolt  impulses, 
construction  was  simplified  because  of  its  smaller  physical  size.  Additionally,  the 
proposed  draft,  lEC  60-2,  permits  the  qualification  of  a larger  divider  by  the 
comparative  measurement  technique  with  a smaller  reference  divider  at  20  % of  the 
rated  voltage  of  the  larger  one  [107].  A reference  divider  rated  at  300  kV  can  therefore 
be  used  to  compare  with  dividers  having  ratings  of  up  to  1500  kV  [107].  The  NIST4 
divider  has  a nominal  resistance  of  10  kQ,  which  is  typical  for  most  resistive  dividers. 
A divider  ratio  of  1000:1  was  selected  rather  than  a larger  one  to  give  a signal-to-noise 
ratio  that  wa^  high  enough  that  degradation  of  the  divider  output  signal  by  pickup 
of  radiated  noise  would  be  minimized.  Both  digital  recorders  and  digital  storage 
oscilloscopes  require  lower  signal  levels  than  this  ratio  provides  at  voltages  exceeding 
200  kV  and  further  attenuation  is  needed.  Attenuators  are  used  at  the  recorder  end 
of  the  signal  cables.  The  divider  and  attenuators  must  have  adequate  bandwidth 
for  the  frequency  range  of  interest.  A resistive  divider  must  have  sufficiently  fast 
transient  response  to  precisely  measure  steep  impulse  voltages  such  as  the  1.2  /zs  front 
time  standard  lightning  impulse.  It  must  have  a stable  divider  ratio  for  the  voltage, 
temperature  and  humidity  range  of  use.  The  voltage  coefficient  and  environmental 
effects  must  be  negligible.  Mechanical  effects  that  unfortunately  have  been  ignored 
in  many  cases  are  also  important.  The  divider  ratio  must  not  be  affected  by  these 
factors  as  well. 
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Figure  50.  Design  of  NIST4. 


4.2.2  Configuration  and  Construction 

The  dimensions  of  the  NIST4  divider  are  shown  in  figure  50.  The  high  voltage 
arm  consists  of  a noninductive  wire-wound  resistor  immersed  in  an  oil-filled  column, 
permitting  a compact  configuration.  The  high-voltage  arm  can  be  physically  shorter 
in  oil  than  in  air  since  the  flashover  voltage  along  the  arm  when  immersed  in 
transformer  oil  is  much  higher  than  in  air.  The  high-voltage  arm  is  only  0.3  m in 
length  while  those  of  most  conventional  300  kV  free-standing  dividers  are  typically 
more  than  three  times  that  length.  The  diameter  of  the  high-voltage  arm  is  only 
2.54  cm,  which  is  also  smaller  than  most  other  dividers  having  comparable  voltage 
ratings. 

There  are  several  advantages  to  make  the  divider  so  compact.  The  compact 
configuration  makes  the  high-voltage  arm’s  stray  capacitance  and  self-inductance 
much  less  than  what  is  possible  when  air  is  used  as  the  insulating  medium.  Lower 
stray  capacitance  gives  the  divider  a faster  transient  response.  Because  the  high- 
voltage  arm  is  shorter,  it  has  a more  uniform  voltage  distribution  along  its  length 
than  otherwise  possible,  which  provides  a more  linear  scaling  of  the  voltage  at  the 
divider  output. 
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A high  voltage  measuring  system  consists  of  more  than  just  the  divider,  however. 
A high  voltage  lead  is  used  to  connect  the  test  object  to  the  divider,  and  this  lead 
makes  the  step  response  for  the  entire  measurement  system  slower  than  that  for  the 
divider  alone.  To  reduce  eifects  arising  from  other  grounded  or  charged  objects  on 
the  measurements,  the  lEC  Standard  recommends  that  the  lead  be  at  least  as  long 
as  1.5  times  the  height  of  the  divider  [103].  A shorter  divider  means  that  a shorter 
high  voltage  lead  may  be  used  without  introducing  significant  proximity  effects,  and 
the  measurement  system  will  have  a faster  overall  transient  response. 

Other  advantages  of  using  oil  insulation  are  better  heat  dissipation  and  reduction 
of  environmental  effects.  A high-voltage  arm  immersed  in  oil  has  much  better  heat 
conduction  than  one  in  air  so  that  temperature  effects  are  smaller.  Leakage  currents 
on  the  surface  of  the  high-voltage  arm  are  negligible  and  unaffected  by  humidity, 
since  the  arm  is  in  oil.  The  column  containing  the  oil  and  the  high-voltage  arm  is 
placed  on  a grounded  metal  cylinder  so  that  the  high-voltage  arm’s  resistance  and  the 
voltage  divider  ratio  would  be  unaffected  by  column’s  surface  resistance  regardless  of 
how  high  the  humidity  in  the  test  area  is. 

Most  external  mechanical  forces  act  on  the  oil-filled  column  instead  of  acting  directly 
on  the  high-voltage  arm  itself  as  is  the  case  with  most  free-standing  resistive  dividers. 
In  this  way,  effects  of  mechanical  strain  on  the  resistance  of  the  high-voltage  arm  are 
minimized.  It  is  shown  later  in  this  section  how  important  these  effects  are  to  the 
value  of  the  resistance. 

It  is  important  the  there  be  no  significant  corona  from  the  divider  or  leads  when  high 
voltage  is  applied.  It  was  found  that  the  high-voltage  arm  immersed  in  oil  was  corona- 
free  under  300  kV  impulse  by  observation  with  a micro-channel  plate  image  intensifier 
camera  in  a darkened  room.  Different  high-voltage  leads  were  also  investigated  in  this 
manner.  Leads  of  2.54  cm  or  smaller  diameter  were  seen  to  have  substantial  corona 
at  that  voltage.  Although  leads  having  diameters  exceeding  10.2  cm  (4  in)  were  seen 
to  be  corona- free,  it  is  not  convenient  to  use  such  large  leads.  It  was  found  that  a 
1 mm  diameter  lead  in  an  oil- filled  tube  was  also  corona- free.  Use  of  this  type  of  lead 
can  provide  a corona-free  impulse  measuring  system,  at  least  to  the  300  kV  level. 

A specially-designed  ring  is  employed  at  the  bottom  of  the  high-voltage  arm  of  the 
NIST4  divider.  Unlike  similar  rings  used  for  capacitive  grading  with  most  300  kV 
resistive  dividers,  the  ring  used  was  specifically  designed  to  shield  the  lower  part  of 
the  high-voltage  arm  to  reduce  noise  pickup  instead  of  providing  improved  voltage 
distribution  along  the  high-voltage  arm.  It  was  seen  from  calculations  based  on  circuit 
models  for  10  kf)  resistive  dividers  that  high-frequency  noise  on  the  high  voltage  leads 
and  top  rings  of  the  dividers  should  be  greatly  damped  by  the  10  kO  resistance  of 
the  high-voltage  arm.  In  step  response  measurements  with  most  resistive  dividers, 
however,  significant  noise  still  appears.  It  was  found  through  comparison  of  different 
shielding  configurations  that  most  of  the  noise  in  the  measured  signal  is  picked  up  by 
the  lower  part  of  high-voltage  arm  where  the  resistance  may  not  be  high  enough  to 
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Figure  51.  Step  Response  With  and  Without  Bottom  Rings. 

attenuate  radiated  noise.  It  is  therefore  necessary  to  have  an  effective  shielding  ring  at 
the  bottom  of  the  high-voltage  arm  to  eliminate  this  noise  and  provide  an  undistorted 
signal.  Figure  51  shows  the  step  response  with  and  without  several  different  bottom 
rings.  The  final  version  of  the  ring  used  in  the  NIST4  divider  was  located  very  close 
to  the  bottom  part  of  high-voltage  arm  to  obtain  the  best  shielding.  A compromise 
was  made  to  get  both  minimal  signal  distortion  and  fast  response.  A larger  ring  of  a 
small  diameter  surrounding  the  bottom  end  of  the  high-voltage  arm  would  have  been 
more  effective  for  shielding,  but  the  increased  capacitance  between  the  high-voltage 
arm  and  the  ring  would  produce  a slower  transient  response.  Additionally,  care  must 
be  taken  to  avoid  flashover  between  the  ring  and  high-voltage  arm.  The  ring  used 
here  was  small,  having  an  inside  diameter  of  5.1  cm  (2  in)  and  a height  of  5.1  cm  (2  in). 
The  NIST4  divider’s  step  response  shown  in  figure  52  is  faster  and  less  distorted  by 
noise  than  those  of  the  other  300  kV  free-standing  dividers,  such  as  the  NIST3  divider 
which  has  a large  big  bottom  ring  of  1 meter  in  diameter. 

There  are  two  parts  that  shield  the  top  of  the  NIST4  divider.  One  is  a sphere  10.2  cm 
(4  in)  in  diameter  on  top  of  the  divider  that  is  used  to  reduce  electrical  field  stresses 
and  eliminate  corona.  The  second  shield  is  a copper  cap  located  inside  the  oil  column 
to  provide  the  fast  response  and  better  voltage  distribution  along  the  high-voltage 
arm.  The  cap  was  empirically  designed  so  that  the  fastest  step  response  possible 
was  achieved  with  acceptable  overshoot.  The  edge  of  the  cap  was  made  of  a 1.27  cm 
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Figure  52.  NIST4  Step  Response. 
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(0.5  in)  diameter  copper  tube  to  ensure  that  corona  would  not  be  produced  at  the 
rated  voltage. 

The  resistive  wire  used  in  the  high-voltage  arm  of  the  NIST4  divider  was  an  alloy-type 
wire  having  a temperature  coefficient  of  -3  ppm  per  degree  centigrade.  Its  diameter  of 
0.11mm  (0.0045  in)  is  thicker  than  that  used  in  most  other  resistive  dividers,  giving 
the  high-voltage  arm  a greater  heat  capacity  and  lower  temperature  increase  after 
application  of  high  voltage  impulses.  The  low  temperature  coefficient,  larger  diameter 
of  the  wire,  and  the  thermal  conductivity  of  the  oil  keep  the  divider’s  resistance  very 
stable  during  environmental  temperature  changes  or  when  several  applications  of  full 
lightning  impulses  are  made.  The  resistance  of  the  high-voltage  arm  was  measured 
with  a resistance  bridge  to  be  10.074  kH  at  room  temperature.  Immediately  following 
the  application  of  ten  consecutive  standard  lightning  impulses  of  325  kV  peak  voltage 
over  a period  of  75  seconds,  the  measured  resistance  was  10.075  kH.  The  difference 
in  resistance  can  be  considered  a negligible  change.  The  low-voltage  arm’s  measured 
resistance  had  the  same  value  of  12.513  0 both  before  and  after  the  ten  impulses. 

The  temperature  coefficient  of  resistance  for  the  high-voltage  arm  depends  not  only 
upon  heating  of  the  wire  but  also  upon  the  changing  mechanical  stresses  of  the 
wire.  An  interesting  experiment  investigating  temperature  and  mechanical  effects  was 
conducted.  Three  sample  resistors  of  the  same  diameter  were  made  with  resistive  wire 
and  placed  into  an  oven.  The  first  sample  was  a simple  wire  loop.  The  second  was 
wound  on  a polyethylene  rod  and  the  third  was  wound  on  a machinable  glass  ceramic 
rod.  Table  1 shows  the  change  in  resistance  after  heating  the  resistors  from  20°  C to 
60° C.  The  resistance  of  the  sample  wound  on  the  polyethylene  rod  had  been  changed 
0.08%  after  heating  to  41°C.  It  is  believed  that  this  is  not  due  to  the  temperature 
coefficient  of  the  winding,  but  rather  to  the  thermal  expansion  of  the  rod  during 
heating  which  resulted  in  a strain-induced  resistance  change.  Changes  in  resistance  of 
this  magnitude  are  unacceptable  for  reference  measurement  systems.  The  machinable 
glass  ceramic  had  a very  good  resistance  stability  at  high  temperature.  Table  2 shows 
the  linear  thermal  coefficient  of  expansion  for  a few  materials  that  are  candidates  for 
substrates  on  which  to  wind  the  resistive  wire.  The  glass  ceramic  has  an  expansion 
coefficient  ten  times  smaller  than  the  other  three  commonly-used  materials.  As  a 
result,  it  was  selected  for  NIST4  divider  as  a base  for  the  high-voltage  arm. 

The  high-voltage  arm  was  wound  in  the  typical  manner  for  resistive  dividers.  Two 
layers  of  windings  are  counterwound  to  minimize  inductance.  The  inner  one  is  wound 
in  clockwise  fashion  and  the  outer  one  in  counterclockwise  fashion.  The  difference 
between  the  windings  of  the  NIST4  divider  and  most  air-insulated  dividers  is  that 
the  rod  on  which  the  wire  is  wound  has  a much  smaller  diameter  and  separation 
between  turns.  There  are  65.7  turns  per  centimeter  (167  turns  per  inch).  From 
prior  experience  [110],  it  is  known  that  the  smaller  the  rod  diameter  and  the  more 
closer  the  spacing  of  the  windings  are,  the  smaller  the  net  magnetic  flux  produced 
by  the  two  windings  will  be.  These  are  important  considerations  in  minimizing  the 
residual  inductance  of  the  wound  resistor.  The  windings  begin  on  a brass  end  cap  and 
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Table  1.  Resistance  Change  of  Resistor  Samples  During  Heating. 


T(°C) 

Simple 

Wire  Loop(A;n) 

Polyethylene 

Rod(n) 

Glass  Ceramic 
Rod{kQ.) 

20 

0.59433 

2.5204 

10.0777 

41 

- 

2.5224 

- 

62 

0.59431 

- 

10.0784 

Table  2.  Linear  Thermal  Coefficients  of  Expansion  of  Candidate  Substrates. 
Polyethylene  Polytetrafluoroethylene  Nylon  Glass  Ceramic 

1.7  X 10-^  1.0  X lO-'^  1.5  X 10-^  9.4  x 10"® 


continue  along  the  ceramic  rod  to  the  opposite  brass  end  cap.  After  the  first  layer 
was  wound,  an  epoxy  coating  was  applied  to  provide  insulation  between  the  windings. 
The  second  layer  of  wire  was  then  counterwound  on  top  of  the  first.  The  epoxy  layer 
used  between  the  windings  is  thin,  compatible  with  transformer  oil,  and  long-lasting 
at  temperatures  below  105° C.  Since  it  is  thin,  the  diameters  of  the  inner  and  outer 
windings  are  nearly  the  same,  which  also  helps  to  reduce  the  overall  inductance  of 
the  resistor.  After  both  windings  were  complete,  their  ends  were  secured  with  screws 
to  the  brass  end  caps  to  ensure  that  there  is  no  mechanical  stress  on  the  wires  caused 
by  moving  the  divider. 

The  low- voltage  arm  of  the  NIST4  divider  has  low  inductance,  high  wattage,  and  good 
shielding.  It  comprises  sixteen  20017,  IW  tin  oxide  resistors  in  parallel.  There  are  two 
shields  for  the  low- voltage  arm,  an  inner  cylindrical  copper  shield  and  an  outer  steel 
cylinder.  Together  they  provide  excellent  shielding  to  electromagnetic  interference. 
The  measured  step  response  for  the  low- voltage  arm  is  shown  in  figure  53,  from  which 
the  risetime  is  calculated  to  be  less  than  1 ns. 


4.3  The  Unit  Step  Response  (USR)  and  Its  Application 

4.3.1  Step  Response  Measurement 

Step  responses  of  the  NIST4  divider  were  taken  using  different  damping  resistors  in 
the  circuit.  The  geometry  used  for  the  step  response  measurements  is  the  horizontal 
loop  configuration  shown  schematically  in  figure  54.  In  this  configuration,  the 
mercury-wetted  relay  used  to  generate  the  voltage  step  is  mounted  on  a support 
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Figure  53.  Step  Response  of  the  Low- Voltage  Side  of  NIST4 
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Figure  54.  Horizontal  Loop  Configuration  for  Step  Response  Measurement 

at  a height  of  about  1 m above  the  ground  plane  and  connected  to  ground  with  a 
wide  copper  strip  labeled  G in  the  figure.  As  mentioned  previously,  the  step  response 
is  changed  by  the  diameter  and  length  of  the  leads  connecting  the  step  generator  to 
the  divider.  A comparison  of  the  step  responses  for  several  different  leads  is  shown 
in  figure  55.  As  seen  from  the  figure,  the  thinner  and  longer  the  lead,  the  worse  the 
step  response,  having  longer  risetimes,  greater  overshoot,  and  more  oscillations. 

Table  3 compares  the  step  response  parameters  T„  (experimental  response  time), 
(partial  response  time),  and  To  (initial  distortion  time)  for  the  NIST4,  NIST3,  and 
NISTl  dividers,  taken  with  different  damping  resistances,  Rd. 


step  Response 


4.3  The  Unit  Step  Response  (USR)  and  Rs  Application 


83 


Figure  55.  Comparison  of  Step  Responses  for  Different  Lead  Lengths 
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Figure  56.  Step  Responses  of  NIST4,  NIST3,  and  NISTl 

The  NIST3  and  NISTl  dividers  are  conventional  air-insulated,  free-standing  resistive 
devices  rated  for  300  kV  and  1 MV  impulses,  respectively.  Figure  56  shows  a 
comparison  of  the  step  responses  for  these  dividers.  Compared  with  other  300  kV 
dividers  used  in  reference  systems  such  as  NIST3,  the  NIST4  divider’s  step  response 
was  much  faster  and  less  distorted.  For  the  undamped  cases,  (i.e.,  no  external 
damping  resistor  in  the  step  response  circuit),  the  10%-90%  risetime  of  the  step 
responses  for  the  three  dividers  is  summarized  in  table  4.  The  NIST4  divider  clearly 
has  a much-improved  risetime. 


4 >3. 2 Convolution  Calculation 

The  convolution  technique  can  be  used  with  the  measured  step  response  to 
predict  the  output  of  an  impulse  measurement  system  for  an  assumed  analytical 
or  numerical  input.  This  approach,  performed  in  the  time  domain,  has  been 
presented  previously  [111].  A frequency-domain  convolution  method  for  estimating 
measurement  errors  is  now  described. 

For  a linear  system  with  an  input  step  voltage  L(t)  and  an  impulse  response  H(t), 
the  step  response  G(t)  is; 


G{t)  = H{t)^L{t) 


(1^) 
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Table  3.  Step  Response  Parameters  for  Three  Reference  Dividers 


Divider 

Rd 

T„ 

Tc 

To 

Overshoot 

(fi) 

(ns) 

(ns) 

(ns) 

(%) 

NIST4 

0 

0.3 

5.3 

0.6 

37.5 

250 

5.5 

8.8 

0.6 

15.2 

500 

14 

14.3 

0.2 

1.7 

NIST3 

0 

19 

17 

2.6 

24 

150 

23 

20 

1.8 

11 

350 

28 

29 

1.3 

3 

550 

36 

36 

0.8 

1 

NISTl 

0 

23 

27 

2.9 

39 

9 

150 

33 

36 

2.7 

14 

250 

41 

44 

1.7 

6 

400 

52 

53 

1.8 

2 

Table  4.  Rise  Time  (ns)  Comparison  for  Three  Reference  Dividers 


NIST4 

NIST3 

NISTl 

11 

40 

52 
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Table  5.  Calculated  Peak  Voltage  Outputs  of  NIST4  Divider 


Rd 

m 

Full  lightning 
(front-time 
0.915  fis) 

Chopped  lightning 
(Chopped  at  88% 
of  the  peak,  Time- 
to-chop  0.566  ^s) 

Chopped  lightning 
(Chopped  at  98% 
of  the  peak,  Time- 
to-chop  0.957  fis) 

0 

1.0000 

0.9991 

0.9993 

250 

1.0000 

0.9948 

0.9990 

500 

1.0000 

0.9914 

0.9983 

where  the  represents  the  convolution  integral.  Similarly,  for  an  input  signal  X(t) 
the  corresponding  output  signal  Y(t)  is: 


Y{t)  = H(t)  * X{t) 


(15) 


In  frequency  domain,  if  G(t),  L(t),  H(t),  X(t),  and  Y(t)  have  Fourier  transforms  G(f). 
L(f),  H(f),  X(f),  Y(f)  respectively!  then 


From  (16)  and  (17), 


G{f)  = H{f)  ■ Lif) 

(16) 

Yif)  = H{f)  ■ X(f) 

(!■) 

X{f)  ■ G{f) 

Hf) 

(IS) 

From  this  formula,  if  a system’s  step  response  is  known,  one  can  predict  an  output 
for  any  given  input  signal  by  taking  the  inverse  Fourier  transform  of  equation  (IS). 
Full  standard  lightning  impulses  along  with  the  chopped  lightning  impulses  shown  in 
figure  57  were  used  as  input  signals  to  calculate  corresponding  output  of  the  NIST4 
divider.  Table  5 gives  some  of  the  peak  values  of  the  output  waveforms  for  normalized 
inputs. 


The  data  of  table  5 indicate  that  the  NIST4  divider  can  measure  the  peak  voltage 
of  full  lightning  impulses  with  negligible  error.  The  data  also  show  that  the  NIST4 
divider  has  less  than  1 % error  in  peak  voltage  measurement  of  front-chopped  lightning 
impulses. 
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Figure  57.  Typical  (a)  Full  and  (b)  Chopped  Lightning  Impulses 
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4.4  Conclusion 

A new,  compact  resistive  voltage  divider  has  been  designed,  constructed,  and 
characterized  for  use  as  part  of  an  impulse  voltage  measurement  system.  Although 
rated  for  impulse  voltages  up  to  300  kV,  it  stands  only  O.Sm  tall.  It  has  significant 
advantages  over  other  free-standing  dividers  because  of  its  small  physical  size,  oil 
insulation,  and  shielding.  The  small  physical  size  gives  it  a response  time  that  is 
nearly  a factor  of  four  times  faster  than  references  dividers  with  comparable  voltage 
ratings  while  the  oil  insulation  provides  much  better  insensitivity  to  environmental 
effects,  such  as  temperature  changes  and  humidity  than  air-insulated  diidders  do. 
A specially-designed  shield  for  the  bottom  of  the  high-voltage  side  of  the  divider 
further  reduces  distortion  of  the  divider  output  by  minimizing  the  pickup  of  spurious 
signals.  The  divider’s  measured  step  response  of  10  ns  is  much  faster  than  comparable 
free-standing  dividers.  The  NIST4  divider  represents  significant  improvement  in  the 
measurement  of  high  voltage  standard  lightning  impulses  with  free-standing  resistive 
voltage  dividers. 
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